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ABSTRACT 


In  this  thesis,  a  perfonnance  analysis  and  improvement  of  a  phase  sampling  inter¬ 
ferometer  antenna  system  based  on  the  Robust  Symmetrical  Number  System  (RSNS)  in 
the  presence  of  noise  is  investigated.  Previous  works  have  shown  that  the  RSNS-based 
DF  technique  can  provide  high  bearing  resolution  with  a  minimum  number  of  antenna 
elements.  However,  the  previous  experimental  data  showed  significant  deviation  from  the 
theoretical  results  expected  due  to  imperfections,  errors,  and  noise.  Therefore,  an  additive 
Gaussian  noise  model  of  RSNS-based  DF  was  established  and  simulated.  Simulation  re¬ 
sults  show  that  the  presence  of  noise  distorts  the  signal  amplitudes  used  in  the  RSNS 
processor  and  causes  degradation  of  the  angle-of-arrival  estimates.  A  performance  analy¬ 
sis  was  undertaken  by  first  introducing  the  quadrature  modulation  configuration  into 
RSNS-based  DF  system,  which  provided  a  digital  antenna  approach  for  more  flexibility 
in  the  signal  processing.  With  a  digital  approach,  variable  resolution  signal  preprocessing 
can  be  employed,  using  a  virtual  channel  concept.  The  virtual  channel  concept  changes 
moduli  values  without  changing  the  actual  physical  antenna  element  spacing.  This  attrac¬ 
tive  property  allows  the  RSNS  algorithm  to  be  implemented  into  existing  antenna  arrays 
and  only  requires  modifying  the  antenna  signal  processor.  Computer  simulation  results 
showed  that  the  proposed  method  can  successfully  improve  the  system  perfonnance  and 
also  mitigate  the  effects  of  noise. 
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EXECUTIVE  SUMMARY 


Radio  direction  finding  (RDF)  systems  acquire  the  angle  of  arrival  (AO A)  of  an 
incident  electromagnetic  wave  (EM)  by  extracting  the  phase  or  time  of  arrival  informa¬ 
tion  from  the  intercepted  waveform.  The  AOA  data  is  prerequisite  information  in  many 
modem  technical  applications,  especially  in  crisis  reaction  areas.  One  of  the  more  signifi¬ 
cant  examples  in  the  military  sector  is  Early  Warning  Threat  Detection,  [1]  which  refers 
to  methods  that  detect  and  identify  the  presence  of  an  imminent  threat  of  attack  and  im¬ 
mediately  take  countermeasures.  The  design  of  a  RDF  device  may  be  divided  into  three 
categories:  amplitude  comparison,  phase  delay,  and  time  delay.  The  phase  delay  tech¬ 
nique  usually  can  generate  an  AOA  resolution  of  ±  1  degree,  which  satisfies  modem  elec¬ 
tronic  warfare  (EW)  requirements.  Thus,  most  radio  direction  finding  systems  are  based 
on  phase-delay  techniques. 

The  direction  finder  based  on  the  phase  delay  technique  increases  element  spacing 
to  improve  the  system’s  resolution;  however,  it  poses  an  ambiguity  problem.  It  exists 
when  a  direction  finder  is  composed  of  only  two  antenna  elements  whose  separation  is 
greater  than  half  of  the  EM  wavelength.  There  exists  a  variety  of  schemes  to  overcome 
the  problem,  such  as  The  Maximum  Likelihood  (ML)  Method,  Eigenspace  Methods,  and 
Estimate  of  Signal  Parameters  via  Rotational  Invariance  Technique  (ESPRIT).  However, 
due  to  the  corresponding  computational  complexity,  few  practical  systems  have  adopted 
these  high  resolution  techniques. 

This  thesis  studied  a  new  phase-sampled  linear  interferometer  technique  that  can 
be  easily  implemented  to  provide  a  high-resolution,  small  baseline  array.  The  approach 
relies  on  the  unique  properties  of  the  robust  symmetrical  number  system  (RSNS)  [7-10], 
which  decomposes  the  spatial  filtering  operation  into  a  number  of  parallel  sub-operations 
(moduli)  that  are  of  small  computational  complexity.  Each  sub-operation  is  considered  as 
a  pre-configured  interferometer  that  symmetrically  folds  the  detected  phase  difference 
with  folding  period  equal  to  2 Nmt  where  i  =  1,2,..., N,  with  N  the  number  of  interfer- 
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ometers  that  are  used  within  the  linear  array,  and  mi  is  the  corresponding  modulus  of  in¬ 
terferometer  i. 

Previous  experimental  results  from  RSNS-based  direction  finders  have  shown  the 
basic  features  of  this  signal  processing  technique.  Since  the  previous  designs  were  under 
noise-free  theoretical  derivations,  the  performances  were  marginal  with  some  large  AOA 
reporting  errors.  Therefore,  the  existing  equations  were  reexamined  and  reconstructed 
from  the  standpoint  of  an  additive  Gaussian  noise  model.  The  system  performance  as¬ 
sessment  for  variable  SNR  was  investigated  and  some  fundamental  relationships  between 
SNR  and  the  AOA  resolution  are  presented. 

This  study  focused  on  developing  a  preprocessing  step  between  the  elements  and 
the  original  RSNS  processing  circuit  to  improve  the  perfonnance  of  the  previous  RSNS- 
based  direction  finder.  Periodic  element  spacing  in  RSNS  algorithm  concept  was  investi¬ 
gated,  which  has  the  benefit  to  be  implemented  into  existing  antenna  arrays  with  simply 
modifying  the  antenna  signal  processor. 

Finally,  a  variable  resolution  signal  processing  approach  based  on  the  virtual 
channel  concept  was  introduced  that  provides  an  increased  resolution  when  the  SNR  is 
high.  Simulation  results  corroborated  the  analysis  and  clearly  demonstrated  the  signifi¬ 
cant  reduction  in  AOA  estimation  error  that  can  be  achieved  with  this  new  design  con¬ 
cept. 
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I.  INTRODUCTION 


A.  RADIO  DIRECTION  FINDING  SYSTEM 

Radio  direction  finding  (RDF)  systems  acquire  the  angle  of  arrival  (AO A)  of  an 
incident  electromagnetic  (EM)  wave  by  extracting  the  phase  or  time  of  arrival  informa¬ 
tion  from  the  intercepted  waveform.  The  AOA  data  is  prerequisite  information  in  many 
modem  technical  applications,  especially  in  crisis  reaction  areas.  One  of  the  more  signifi¬ 
cant  examples  in  the  military  sector  is  Early  Warning  Threat  Detection,  [1]  which  refers 
to  methods  that  detect  and  identify  the  presence  of  an  imminent  threat  of  attack  and  im¬ 
mediately  take  countermeasures.  Therefore,  numerous  techniques  have  been  developed  to 
satisfy  a  wide  variety  of  requirements. 

The  design  of  a  RDF  device  may  be  divided  into  three  categories,  depending  on 
the  fundamental  measurement  principles  [2]: 

Category  I:  Amplitude  comparison. 

Category  II:  Phase  delay. 

Category  III:  Time  delay. 

A  Category  I  system  converts  the  amplitude  response  of  the  antenna  elements  into 
voltages  in  order  to  obtain  the  AOA  information.  Figure  1  depicts  the  two  types  of  an¬ 
tenna  amplitude  techniques  used  to  obtain  AOA  information.  They  are  referred  to  as  the 
direct  and  comparative  amplitude  techniques.  Here,  9  represents  the  AOA  to  be  meas¬ 
ured. 

For  Categories  II  and  III,  at  least  two  separate  antenna  elements  with  predeter¬ 
mined  space  between  them  are  required.  Figure  2  depicts  the  basic  configuration  for  the 
phase-delay  and  time-delay  techniques.  We  will  assume  a  single  frequency  time- 
harmonic  waveform.  Since  the  electromagnetic  wave  experiences  different  propagation 
distances  from  an  emitter  to  the  individual  antenna  elements  (except  in  the  broadside 
case),  it  induces  different  phase  terms  at  the  antenna  terminals.  Equivalently,  the  wave 
arrives  at  each  antenna  element  at  a  different  time.  The  phase  difference  or  time  differ- 
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ence  is  dependent  on  the  element  spacing,  wavelength,  and  AOA.  Since  the  first  two  fac¬ 
tors  are  known,  the  system  can  compute  the  AOA  information. 


Figure  1. 


Horizontal  E-field  Dipole 

(a)  (b) 

(a)  Direct  and  (b)  Comparative  Amplitude  Techniques.  (After  [2].) 


Figure  2.  Basic  Phase-Delay/Time-Delay  Techniques.  (After  [2].) 

Phase  comparison  systems  usually  can  generate  an  AOA  resolution  of  ±  1  degree, 
which  satisfies  modem  electronic  warfare  (EW)  requirements  [1].  Thus,  most  radio  direc¬ 
tion  finding  systems  are  based  on  phase-delay  techniques.  However,  to  improve  the  sys¬ 
tem’s  resolution,  the  element  spacing  must  be  increased,  which  poses  an  inherent  prob¬ 
lems  of  direction  ambiguities.  They  exist  for  the  phase-delay  techniques  when  an  inter¬ 
ferometer  is  composed  of  only  two  antenna  elements  whose  separation  is  greater  than 

half  of  the  signal’s  wavelength.  These  ambiguities  cause  error  estimates  in  the  angle  of 
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arrival.  There  exists  a  variety  of  schemes  to  resolve  the  ambiguities,  several  of  which  are 
reviewed  in  the  following  paragraphs. 

The  Maximum  Likelihood  (ML)  Method  [3]  estimates  the  AOA  from  a  given  set 
of  sampled  phased  data  by  maximizing  the  log-likelihood  function.  The  likelihood  func¬ 
tion  is  the  joint  probability  density  function  of  the  sampled  data,  given  the  AOA,  and 
viewed  as  a  function  of  the  desired  variables  (the  AOA  in  this  case).  The  ML  methods 
yield  nearly-optimal  performance.  Even  when  the  signal-to-noise  ratio  (SNR)  and  the 
number  of  samples  are  small,  they  still  have  better  performance  compared  to  other  meth¬ 
ods.  ML  methods,  however,  require  significant  computational  resources. 

Eigenspace  Methods  rely  on  a  property  of  the  array  covariance  matrix.  The  space 
spanned  by  its  eigenvectors  may  be  partitioned  in  two  subspaces,  namely,  the  signal  sub¬ 
space  and  the  noise  subspace  when  sources  are  incoherent.  The  steering  vectors  corre¬ 
sponding  to  the  directional  emitters  are  orthogonal  to  the  noise  subspace.  In  principle  the 
Eigenspace  methods  search  for  directions  such  that  the  steering  vectors  associated  with 
these  directions  are  orthogonal  to  the  noise  subspace  and  are  contained  in  the  signal  sub¬ 
space.  The  Multiple  Signal  Classification  (MUSIC)  [4]  algorithm  is  an  eigenspace 
method.  Many  improvements  have  been  proposed  to  achieve  better  estimation  perform¬ 
ance,  such  as  Root-MUSIC  [5].  The  Root-MUSIC  algorithm  solves  for  the  root  of  a 
polynomial  instead  of  perfonning  the  conventional  MUSIC  search,  and  exhibits  better 
performance  and  a  more  reasonable  computational  load. 

A  third  method,  Estimate  of  Signal  Parameters  via  Rotational  Invariance  Tech¬ 
nique  (ESPRIT)  [6]  has  a  lower  computation  and  storage  requirement  than  other  eigen¬ 
space  techniques  such  as  MUSIC.  This  advantage  is  achieved  by  using  two  identical  ar¬ 
rays,  in  the  sense  that  array  elements  need  to  form  matched  pairs  with  identical  displace¬ 
ment  by  the  same  distance  and  in  the  same  direction  relative  to  the  first  element.  An  array 
geometry  must  be  designed  to  have  this  property.  Due  to  the  corresponding  computa¬ 
tional  complexity,  few  practical  systems  have  adopted  these  high  resolution  techniques. 

This  research  studied  a  new  phase  sampled  linear  interferometer  technique  that 
can  be  easily  implemented  to  provide  a  high-resolution,  small  baseline  array.  The  ap¬ 
proach  relies  on  the  unique  properties  of  the  robust  symmetrical  number  systems  (RSNS), 
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[7-10]  which  decompose  the  spatial  filtering  operation  into  a  number  of  parallel  sub¬ 
operations  (moduli)  that  are  of  small  computational  complexity.  Each  sub-operation  is 
considered  as  a  pre-configured  interferometer  that  symmetrically  folds  the  detected  phase 
difference  with  folding  period  equal  to  2 Nmi  where  i  =  1,2,..., TV,  with  N  the  number  of 

interferometers  that  are  used  within  the  linear  array  and  mi  is  the  corresponding  moduli 
of  interferometer  i. 

The  RSNS  technique  can  be  applied  to  either  the  amplitude  or  phase  at  the  ele¬ 
ment  output.  The  parallel  use  of  phase  waveforms  increases  the  antenna  resolution  with¬ 
out  increasing  the  folding  rate  of  the  system.  Small  amounts  of  analog-to-digital  process¬ 
ing  are  required  to  quantize  each  folded  phase  response  so  that  each  sub-operation  only 
requires  a  precision  in  accordance  with  its  modulus.  This  feature  permits  the  system  to 
overcome  phase  errors  induced  by  the  array  environment,  including  mutual  coupling.  A 
much  higher  AOA  resolution  is  achieved  by  using  more  moduli  in  a  system.  Previous  ex¬ 
perimental  results  from  RSNS-based  direction  finder  [7-10]  have  shown  the  effectiveness 
of  this  signal  processing  technique.  The  previous  research  presented  a  “mixed  signal” 
processing  approach.  This  thesis  continued  an  examination  of  the  perfonnance  of  this 
novel  signal  processing  under  noise  environments,  and  how  a  digital  RSNS  architecture 
might  be  implemented. 

B.  PRINCIPAL  CONTRIBUTIONS 

This  thesis  research  focused  on  developing  a  preprocessing  step  between  the  ele¬ 
ments  and  the  original  RSNS  processing  circuit  to  improve  the  perfonnance  of  the  previ¬ 
ous  RSNS-based  direction  finder.  A  new  algorithm  was  established  and  new  features 
were  introduced,  including  the  possibility  of  obtaining  variable  resolution  with  a  fixed 
element  spacing. 

The  previous  RSNS  direction  finders  were  designed  under  noise-free  theoretical 
derivations  and,  therefore,  the  performances  were  marginal  with  some  large  AOA  report¬ 
ing  errors  [7-10].  The  existing  equations  were  re-examined  and  re-constructed  from  the 
standpoint  of  an  additive  Gaussian  (thermal)  noise  model.  Computer  simulations  with 
noise  were  conducted  in  order  to  evaluate  the  previous  system  designs.  The  system  per- 
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formance  assessment  for  variable  SNR  was  investigated.  Some  fundamental  relationships 
between  SNR  and  the  AOA  resolution  are  presented. 

During  the  course  of  the  simulations,  it  was  discovered  that  different  moduli  pairs 
under  the  same  SNR  condition  have  different  performances.  Since  the  selection  of  moduli 
affects  the  system  performance,  the  possibility  of  changing  moduli  values  without  chang¬ 
ing  the  actual  physical  antenna  element  spacing  was  investigated. 

The  possibility  of  employing  equal  spacing  and  still  using  a  RSNS  processing  al¬ 
gorithm  was  also  examined.  Periodic  element  spacing  is  a  very  desirable  characteristic.  If 
the  RSNS  algorithm  can  be  applied  to  periodic  elements,  then  the  concept  can  be  imple¬ 
mented  using  any  conventional  array.  It  would  simply  be  a  matter  of  modifying  the  an¬ 
tenna  signal  processor. 

Finally,  a  variable  resolution  signal  processing  approach  based  on  virtual  channel 
concept  is  introduced  that  provides  an  increased  resolution  when  the  SNR  is  high.  Simu¬ 
lation  results  corroborate  the  analysis  and  clearly  demonstrated  the  significant  reduction 
in  AOA  estimation  error  that  can  be  achieved  with  this  new  design  concept. 

C.  THESIS  OUTLINE 

The  organization  of  this  thesis  is  as  follows.  Chapter  II  reviews  the  basic  operat¬ 
ing  principles  of  RDF  systems,  in  general,  and  the  two-element  phase  sampled  interfer¬ 
ometer,  in  particular.  The  equations  governing  the  sampling  of  an  incident  wave  front  are 
presented  and  the  variables  that  affect  the  output  are  discussed. 

Chapter  III  is  concerned  with  the  RSNS  theory.  The  application  of  RSNS  theory 
to  phase  sampled  interferometry  is  described,  and  the  previous  experimental  results  are 
presented. 

Chapter  IV  addresses  the  establishment  of  a  noise  model  and  the  analysis  of  noise 
effects  in  a  RSNS-based  interferometer.  The  quadrature  demodulation  approach  was  in¬ 
vestigated  to  provide  more  flexibility  on  performance  improvement.  The  possibility  of 
applying  the  RSNS  theory  to  uniformly  spaced  arrays  was  also  studied.  A  variable  resolu¬ 
tion  signal  processing  method  was  proposed  that  provides  a  performance  improvement  in 
the  presence  of  noise. 
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Chapter  V  presents  several  case  studies  and  computer  simulations  using  the  equa¬ 
tions  derived  in  Chapter  IV.  Simulation  results  are  discussed  and  analyzed. 

Chapter  VI  gives  a  summary  of  this  study  and  recommendations  for  further  re¬ 
search. 

The  Appendix  presents  the  MATLAB  codes  used  to  design,  simulate,  and  evalu¬ 
ate  the  RSNS-based  DF  system. 
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II.  PHASE  DELAY-BASED  INTERFEROMETRY 


In  this  chapter,  we  review  the  basic  operating  principles  of  RDF  systems,  in  gen¬ 
eral,  and  the  two-element  phase  sampled  interferometer,  in  particular.  The  equations  gov¬ 
erning  the  sampling  of  an  incident  wave  front  are  presented  and  the  variables  that  affect 
the  output  are  discussed. 

A.  THEORY 

At  the  outset,  it  is  important  to  understand  that  under  ideal  circumstances,  the 
electromagnetic  wave  incident  on  DF  antennas  exhibits  a  locally  plane  wave  structure 
with  linear  polarization  when  the  receiver  is  far  enough  from  the  transmitter  and  the  an¬ 
tenna  interactions  with  the  platform  are  controlled.  An  ideal  direction  tinder  has  the  ca¬ 
pability  of  detecting  signals  over  three-dimensional  space.  It  is  achieved  by  combining 
any  two  orthogonal,  one-dimensional  linear  interferometers  in  a  plane.  Therefore,  a  brief 
description  of  the  underlying  theory  of  a  one-dimensional  linear  interferometer  is  helpful 
in  understanding  the  material  in  later  chapters. 

A  simple  two-element  linear  interferometer  based  on  phase  delay  interferometry 
is  shown  in  Figure  3.  Two  identical  antenna  elements  are  separated  by  a  distance  d  and 
the  plane  wave  impinges  at  an  incident  angle  9  measured  from  the  perpendicular  to  the 
baseline  axis  [7].  The  angle  9  takes  on  values  from  -90°  to  90°. 


Figure  3.  Two-Element  Interferometry.  (After  [7].) 
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Referring  to  Figure  3,  the  incident  wave  arrives  first  at  antenna  element  1  and  then  at 
element  2  after  propagating  an  additional  distance  of  Jsin(<9).  The  signals  received  by 
the  two  antenna  elements  have  the  following  form 

F($,f)  =  Fcos  (ajj  +  yfi+fa)  (2.1) 

where  i  represents  the  antenna  index  number,  V  is  signal  amplitude  response,  OJ  is  the 
carrier  frequency  in  radians/second,  (pl  is  the  phase  term  relative  to  the  origin,  and  y/i  is 

the  phase  tenn  caused  by  the  transmission  line.  The  received  signals  are  fed  into  a  phase 
detector  comprised  of  a  mixer  and  a  low-pass  filter  (LPF).  After  the  phase  detector  mix¬ 
ing  process,  the  received  signal  has  the  form 

y[cos(2^+^+^)+cos(a^+a^)]  (2-2) 

where  Z  and  A  correspond  to  the  sum  and  difference  of  the  phase  terms,  respectively. 

After  LPF  filtering,  the  output  signal  from  the  phase  detector  has  a  DC  voltage  of 
the  form 

CW  =  y«»(A|rt^)  (2.3) 

where  Ay/  comes  from  the  transmission  line  length  differences.  The  term  A<p  represents 
the  phase  differential  between  the  two  received  signals,  expressed  as 

A(p  -  (px- (j)2  =Msin(<9)  (2.4) 

where  k  is  the  wavenumber 


and  X  is  the  wavelength  at  carrier  frequency  (Oc. 

It  is  reasonable  to  assume  that  Ay/  in  Equation  (2.3)  can  be  measured  and  then 
eliminated.  Let  the  tenn  V2  j 2  be  nonnalized  to  unit  amplitude,  so  that  we  obtain  the  fol¬ 
lowing  final  simple  output  signal  from  the  mixing  process: 
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Vout  ( 0 )  =  cos(A^)  =  cos (kd sin(#)). 


(2.6) 


B.  AMBIGUITY  AND  ARTICULATION 

From  Equation  (2.4),  the  phase  differential,  A (j),  repeats  itself  with  a  period  2k. 

Thus  the  phase  measurement  becomes  ambiguous  at  high  frequencies  or  large  spacings. 
Therefore,  a  single  phase  differential  value  corresponds  to  many  angles  of  arrival.  Ambi¬ 
guities  are  present  when  the  distance  between  antenna  elements  is  greater  than  half  of  the 
frequency  wavelength,  or  d  >  X/ 2. 

We  take  two  representative  numerical  examples  to  demonstrate  this  criterion. 

First,  consider  the  case  of  d  -  A/ 2,  which  gives  A0  =  ;rsin(<9).  The  phase  differential 

versus  angle  is  shown  in  Figure  4.  There  is  a  one-to-one  mapping  between  AOA  and 
phase  difference.  There  are  no  AOA  ambiguities  created  by  spacing. 


Figure  4.  Phase  Difference  versus  AOA  for  d  -  A/2.  (After  [8].) 

Next,  consider  the  case  of  d  =  A,  which  gives  A(j)  -  2;rsin(0).  The  phase  differ¬ 
ence  versus  AOA,  shown  in  Figure  5,  suffers  from  ambiguities.  The  circles  shown  in 
Figure  5  have  the  same  phase  differential,  which  leads  to  the  ambiguous  estimation  of 
angle  of  arrival. 
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Figure  5.  Phase  Difference  versus  AOA  for  d  -  A . 


We  can  define  the  AOA  measurement  precision  as  the  rate  of  change  of  the  phase 
differential  versus  AOA,  or  3  (A (/))/d6.  The  unit  is  degrees  per  degree  [2],  Taking  the 
derivative  of  both  sides  of  Equation  (2.4)  yields 


(2.7) 


which  shows  that  the  measurement  precision  is  proportional  to  the  ratio  of  d/X  for  con¬ 
stant  carrier  frequency  and  incident  angle.  Although  a  long  baseline  suffers  from  phase 
ambiguities,  it  increases  the  resolution  of  the  AOA  estimation. 

In  practice,  an  interferometer  DF  system  consists  of  dual  baseline  structures,  so 
that  the  shorter  baseline  is  used  to  overcome  the  ambiguities  and  the  longer  baseline  is 
used  to  achieve  high  resolution  [1]. 

C.  FOLDING  WAVEFORMS 

An  examination  of  Equation  (2.6)  shows  that  the  output  voltage  contains  the  plane 
wave  AOA  information  in  the  fonn  of  a  symmetrical  folding  waveform.  This  feature  is 
illustrated  in  Figure  6  for  the  two  representative  spacing  cases  which  have  been  dis¬ 
cussed.  The  number  of  folding  periods  n  within  the  visible  region  [-7r/2<0<7r/2)  is 
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n  =  d/{A/2)  -  2 d/ A. 


(2.8) 


Angle  of  Arrival  (degrees) 


Figure  6.  Normalized  Output  versus  Incident  Angle  for  (a)  d  -  A/2  (b)  d  =  A.  (After 

mo 

By  relating  the  element  distances  by  a  successive  factor  of  two  as  shown  in  Figure 
7  and,  appropriately  shifting  the  output  folding  waveform  of  each  baseline,  it  is  possible 
to  achieve  a  binary-resolution  interferometer  system  [1],  For  example,  a  4-bit  linear  base¬ 
line  interferometer  consists  of  5  antenna  elements  with  spacings  of  2-1  A,  2°  A,2X  A,  and 
22A.  This  kind  of  harmonic  binary-related  scheme  is  limited  in  its  application  due  to  re¬ 
quirements  of  a  long  baseline  spacing  to  achieve  high  resolution. 

Anti  Ant  2  Ant  3  Ant  4  AntS 


-1/2- 


-2k- 


-41- 


Figure  7.  Four-Bit  Linear  Baseline  Interferometer.  (After  [1].) 

In  summary,  this  chapter  has  discussed  some  basic  direction  finding  techniques 
and  their  limitations.  A  signal  preprocessing  method  employing  the  robust  symmetrical 
number  system  is  introduced  in  the  next  chapter.  It  can  achieve  a  high  resolution  without 
large  element  spacing. 
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III.  ROBUST  SYMMETRICAL  NUMBER  SYSTEM-BASED 

DIRECTION  FINDING 


Since  it  is  desirable  to  obtain  high  resolution  AOA  estimates  without  a  long  base¬ 
line,  a  signal  preprocessing  method  employing  the  robust  symmetrical  number  system 
(RSNS)  was  introduced  in  [7-10].  In  this  chapter,  the  RSNS  theory  is  reviewed  and  pre¬ 
vious  experimental  data  are  presented  and  discussed. 


A.  THE  ROBUST  SYMMETRICAL  NUMBER  SYSTEM 

The  concept  of  RSNS-based  preprocessing  is  that  a  robust  symmetrical  number 
sequence  with  proper  shift  value  is  used  to  represent  a  symmetrical  folding  waveform  out 
of  the  antenna.  Consider  an  A-channel  RSNS,  where  N  >  1 .  The  basic  RSNS  sequence  has 
the  following  form: 


Xm>  =  [0^.  .  dA  C  -  A , ,  m, , . . . ,  m, , AJ] 

N  N  N 


(3.1) 


where  Xm  represents  the  i"'  channel  sequence  (1  <  i  <  N)  with  a  corresponding  modulus 
mr  Equation  (3.1)  originates  from  repeating  each  value  N times  in  the  following  se¬ 
quence: 

[0,  l,2,...,m  — l,m,m  — 1,...,2,1]  (3.2) 


where  m  is  a  positive  integer  (m  >  1). 


It  has  been  shown  that  the  period  of  a  channel  with  modulus  mi  is  [9] 


P  —  9  A Iyyi 

1  RSNS 


(3.3) 


and  the  fundamental  period  of  an  A-channel  RSNS  with  modulus  set  { mx,m2,...,mN }  is 


RRrsns  '2.N\ml,m2,...,mN] 

where  \mvm2,...,mN\  represents  the  least  common  multiple  of  mljn2,...jnN. 


(3.4) 
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The  discrete  states  of  each  channel  can  expressed  as  [9] 


l~s, 

N  _ 

2  Nmi  +N  —  l  +  sj  —  1 
N 


st  <  I  <  Nmi  +  y  + 1 

Nmi  +  y  +  2  <  1  <  2 Nmi  +  y  - 1 


(3.5) 


where  g  is  the  Ith  state  value  of  channel  i,  y  is  the  corresponding  sequence  shift  defined 
as 


y  =  0, 1, 2, •  •  • , iV  —  l(mod  N)  (3.6) 

and  [xj  indicates  the  greatest  integer  less  than  or  equal  to  x.  The  sequence  shift  set 
{y,^,...,^}  must  be  a  complete  residue  system  modulo  N  in  order  to  have  Gray-code 
properties. 

Table  1  lists  discrete  states  within  the  RSNS  for  m,  =  3 (.s',  =  0)  and  m2  =  4 
(s2  =  l).  The  column  vectors,  which  consist  of  the  values  within  each  modulus,  change 
one  at  a  time  between  code  positions  (Gray-code  property). 
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Table  1.  Two-Channel  RSNS  for  Modulus  3  and  4. 


It  is  desired  to  find  the  largest  continuous  set  of  column  vectors  shown  in  Table  1 
that  are  distinct.  This  sequence  of  distinct  vectors  forms  the  unambiguous  output  of  the 
system  and  is  defined  as  the  system  dynamic  range ,  denoted  by  M.  Studies  [11]  have 
shown  that  the  selection  of  the  sequence  shift  (y. )  and  permutations  among  the  moduli 

has  no  effect  on  M,  but  may  produce  different  start  and  stop  points  for  the  vectors  corre¬ 
sponding  to  M.  The  dynamic  range  of  the  2-channel  example  shown  in  Table  1  is 
M  =  15.  That  is,  starting  at  vector  (2,2)  and  ending  at  vector  (3,1),  all  vectors  are 
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unique.  The  symbol  n  re-indexes  the  sequence  of  RSNS  vectors  and  counts  the  first  vec¬ 
tor  in  the  dynamic  range  as  zero  and  the  last  vector  in  the  dynamic  range  as  M  —  1  or  14. 

The  dynamic  range  of  several  particular  moduli  combinations  have  been  formu¬ 
lated  in  closed  form.  For  two  channels,  three  combinations  of  two  relatively  prime 

moduli  m1  and  m2  can  yield  M  as  follows  [10]: 

Case  I:  ml  >  3  and  m2  =  /«,  + 1,  then 

M  =3(ml+m2)-6.  (3.7) 

Case  II:  m{  >  5  and  m2  =  ml  +  2,  then 

M  =  3(mx+m2)-l .  (3.8) 

Case  III:  m1  >  5  and  m2  =  ml  +  C,  and  C  >  3,  then 

M  —  4/«1  +  2m2  —  2.  (3.9) 

Checking  the  example  shown  in  Table  1,  which  fits  Equation  (3.7),  yields  M  =  15. 

So  far,  we  have  reviewed  the  basic  features  of  the  RSNS.  In  the  next  section,  the 
relationship  between  RSNS  theory  and  DF  systems  is  discussed. 


B.  THE  RSNS  ANTENNA  ARCHITECTURE 

There  are  three  steps  to  integrate  RSNS  theory  into  a  direction  finding  antenna 
design;  these  steps  are  presented  in  the  following  paragraphs. 

The  first  is  to  decide  the  physical  antenna  element  spacing  according  to  the  se¬ 
lected  moduli  and  operating  frequency.  Referring  to  Equation  (3.3)  and  the  definition  of 
dynamic  range,  we  define  the  number  of  folds  within  each  modulus  as 


M  _  M 

Prsns  3Nmj 


(3.10) 


Then  combining  Equation  (2.8)  and  (3.10),  we  get  the  corresponding  required  antenna 
spacing  for  modulus  mi 
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(3.11) 


,  X  MX 
a,  =  n,  —  = - . 

2  4  Nmi 

In  practice,  the  antenna  element  pattern  degrades  at  wide  angles  so  that  the  re¬ 
ceived  waveforms  are  not  ideal  over  a  full  180  degrees.  Therefore,  the  RSNS  can  be  re¬ 
mapped  to  the  reliable  portion  of  the  pattern  to  increase  resolution  [7].  Considering  Equa¬ 
tion  (2.6)  with  k  and  Vout  held  constant,  we  have 


d  sin(<9)  =  d'  sin(#') 


(3.12) 


where  0  and  d  correspond  to  the  maximum  field  of  view  (FOV)  and  the  antenna  spacing 
for  the  original  ( unsealed)  configuration.  The  quantities  0'  and  d'  are  the  maximum  FOV 
and  antenna  spacing  for  the  remapped  ( scaled)  configuration.  Rearranging  Equation 
(3.12),  the  scale  factor  is  defined  as 


d'  _  sin(<9) 
d  sin(#') 


(3.13) 


Combining  Equations  (3.11)  and  (3.13),  we  obtain  the  scaled  antenna  spacing 


d'  =  df=h4^f 

1  1  4  m..N 


(3.14) 


Next,  in  an  analog  implementation,  we  need  to  design  the  comparators  which  en¬ 
code  the  normalized  received  waveforms  shown  in  Equation  (2.6)  into  RSNS  form.  The 
number  of  comparators  required  for  each  channel  is  equal  to  the  corresponding  modulus. 
It  has  been  shown  [7]  that  for  a  nonnalized  folding  waveform  amplitude  (-1  <  Vout  <  l), 

the  threshold  for  the  jth  comparator  of  modulus  channel  mi  is 


Vi,mi  =  C0S 


m,-j  + 1/2 


K 


III; 


,  1  <  j  <  m, . 


(3.15) 


Whenever  a  folding  waveform  exceeds  a  comparator  threshold,  the  comparator 
turns  ON.  A  consequence  of  quantizing  a  folding  waveform  into  discrete  states  is  that  the 
RSNS  system  displays  any  signal  that  falls  within  a  discrete  state  bin  as  if  it  arrived  at  the 
bin  center.  Figure  8  illustrates  discrete  states  of  a  two-channel  RSNS  with  symmetrical 
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folding  waveforms  having  moduli  mx  =  3  (5,  =0)  and  m2  =  4(s2  =  l).  The  thresholds 

Vj  m  from  Equation  (3.15),  shown  on  the  vertical  axis,  encode  the  folding  waveforms 

into  the  RSNS.  The  discrete  states  (integers)  shown  above  the  folding  waveforms  repre¬ 
sent  the  number  of  comparators  in  the  ON  state. 

L  M~15  >1  mi=3»V|=0 


0  1  12233221  I  0t»l  122332211  0  0  1  1223  3  22 


Figure  8.  RSNS  Folding  Waveforms  for  mx  =  3 (.s’,  =0)  and  m2  =4 (s2  =  l).  (After 

[8].) 

Finally,  we  have  to  calculate  the  phase  adjustment  tenn  for  each  channel  in  order 
to  map  the  center  of  the  dynamic  range  to  the  response  of  the  antenna  when  the  signal 
AOA  is  at  broadside.  Referring  to  Table  1,  neither  modulus  has  symmetry  about  the  mid¬ 
point  of  the  dynamic  range,  which  is  n  =  7.  Therefore,  a  phase  adjustment  g  is  added  to 
Equation  (2.6)  giving 

VoJ6)  =  cos(kdsm(0)  +  g).  (3.16) 

When  6  =  0,  sin  (6)  =  0,  so  Vout  [6  =  0)  =  cos (g).  The  value  of  g  is  selected  so  that 
Vout  (0-0)  is  encoded  to  the  center  bin  of  the  system  dynamic  range.  In  the  example 
shown  in  Table  1,  Vout  (0  -  0)  should  be  encoded  into  0  for  m]  and  3  for  m2. 

At  this  point,  we  can  obtain  the  relationship  between  the  RSNS  digital  output  and 
the  angle  of  arrival.  Using  this  relation,  we  further  obtain  the  transfer  function  of  the 
processor,  which  is  a  plot  of  true  AOA  versus  estimated  AOA.  Studies  have  shown  that 
the  estimated  AOA  0'u  for  the  u"‘  bin  is  given  as  [7] 
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O'  =  sin 


-1 


f2u  +  \  1A 


MZ,  $ 


(3.17) 


and  the  AO  A  resolution  r  for  the  ulh  bin  is  given  as 


r  =sin 


' lu-M  +  2 ^ 

-sin  1 

2u  —M^ 

M<f  j 

{  M4  ) 

(3.18) 


In  the  next  section,  a  design  example  is  given  to  illustrate  the  important  features 
of  the  RSNS  processing  method. 


C.  DESIGN  EXAMPLE  AND  EXPERIMENTAL  RESULTS 

In  this  section,  a  two-channel,  6-bit  RSNS  antenna  design  originally  given  in  [7] 
is  described  to  illustrate  the  procedures  in  Section  B.  MATLAB  programs  are  also  devel¬ 
oped  in  order  to  provide  essential  computational  modeling  and  simulation.  All  MATLAB 
programs  are  listed  and  explained  in  the  Appendix. 

The  schematic  diagram  shown  in  Figure  9  is  for  the  design  [13]  of  a  6-bit  RSNS- 
based  interferometer  using  m]  =  8  (.s',  =  0)  and  m2  =  1 7  (,s'2  =  1). 


0B 


Figure  9.  Block  Diagram  of  the  RSNS-Based  DF  System.  (From  [13].) 
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The  antenna  elements  are  printed  circuit  dipoles  operating  at  the  frequency  of  8.0 
GHz.  For  these  RSNS  parameters,  we  proceed  to  calculate  the  essential  array  configura¬ 
tion  and  transfer  function  with  the  assistance  of  MATLAB  programs. 

1.  Simulation  Procedure 

Step  1.  By  Equation  (3.9)  we  obtain  the  dynamic  range,  M  =  64.  The  RSNS  vec¬ 
tor  sequence  is  shown  in  Table  2. 


/ 

m2 

n 

/ 

m1 

m2 

n 

/ 

m , 

m2 

n 

21 

6 

10 

0 

43 

5 

13 

22 

65 

0 

2 

44 

22 

5 

10 

1 

44 

6 

13 

23 

66 

1 

2 

45 

23 

5 

11 

2 

45 

6 

12 

24 

67 

1 

1 

46 

24 

4 

11 

3 

46 

7 

12 

25 

68 

2 

1 

47 

25 

4 

12 

4 

47 

7 

11 

26 

69 

2 

0 

48 

26 

3 

12 

5 

48 

8 

11 

27 

70 

3 

0 

49 

27 

3 

13 

6 

49 

8 

10 

28 

71 

3 

1 

50 

28 

2 

13 

7 

50 

7 

10 

29 

72 

4 

1 

51 

29 

2 

14 

8 

51 

7 

9 

30 

73 

4 

2 

52 

30 

1 

14 

9 

52 

6 

9 

31 

74 

5 

2 

53 

31 

1 

15 

10 

53 

6 

8 

32 

75 

5 

3 

54 

32 

0 

15 

11 

54 

5 

8 

33 

76 

6 

3 

55 

33 

0 

16 

12 

55 

5 

7 

34 

77 

6 

4 

56 

34 

1 

16 

13 

56 

4 

7 

35 

78 

7 

4 

57 

35 

1 

17 

14 

57 

4 

6 

36 

79 

7 

5 

58 

36 

2 

17 

15 

58 

3 

6 

37 

80 

8 

5 

59 

37 

2 

16 

16 

59 

3 

5 

38 

81 

8 

6 

60 

3 

16 

17 

60 

2 

5 

39 

82 

7 

6 

61 

3 

15 

18 

61 

2 

4 

40 

83 

7 

7 

62 

4 

15 

19 

62 

1 

4 

41 

84 

6 

7 

63 

mm 

4 

14 

20 

63 

1 

3 

42 

85 

6 

8 

Table  2.  RSNS  Vector  Sequence  Table  for  ml  =  8  and  m2  =  17. (After  [7].) 


Step  2.  With  a  scale  factor  of  £  =  2/^3  ,  by  Equation  (3.13),  the  FOV  changes 
from  ±90°  to  ±60°.  The  corresponding  element  spacings  by  Equation  (3.14)  are 
dl  =  4.3 1  cm  and  d2  =  2.03  cm. 

Step  3.  By  Equation  (3.15)  the  corresponding  thresholds  of  each  channel’s  com¬ 
parators  are  shown  in  Table  3. 
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Comparator 

Channel  ml 

Channel  m2 

Comparator 

Channel  mx 

Channel  m2 

1 

cos  (l5/r/16) 

cos  (33^/34) 

10 

cos  (l5^/34) 

2 

cos  (l3;r/16) 

cos  (31^/34) 

11 

cos  (13^/34) 

3 

cos  (l  lzr/16) 

cos  (29^r/34) 

12 

cos  (l  lzr/34) 

4 

cos  (9;r/16) 

cos  (27^r/34) 

13 

cos  (9^r/34) 

5 

cos  (7/r/16) 

cos  (25^/34) 

14 

COS  (7^/34) 

6 

cos  (5/r/16) 

cos (23^/34) 

15 

COS  (5^/34) 

7 

cos  (3/r/16) 

COS  (2  l/r/34) 

16 

COS  (3^/34) 

8 

cos  (l^r/1 6) 

cos  (l9^r/34) 

17 

cos  (lzr/34) 

9 

cos  (17^/34) 

Table  3.  Comparator  Threshold  Table.  (After  [7].) 


Step  4.  A  phase  adjuster  g.  for  each  channel  is  calculated  using  Equation  (3.16). 
The  normalized  voltage  output  for  channel  m1  after  adding  a  phase  adjustment  is  given  as 

Voun{0)  =  cos{kdlsm(0)  +  7r/4)  (3.19) 


and  for  channel  m2 


Ku, 2  (#)  =  cos  (kd2  sm(0)  +  a/2). 


(3.20) 


Figure  10  depicts  the  folding  output  waveforms  from  Equations  (3.19)  and  (3.20). 
The  corresponding  quantized  folding  output  waveforms  are  plotted  in  Figure  11. 


Figure  10.  Folding  Output  Waveforms  with  Phase  Adjustment. 
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Figure  11. 


Phase  Adjusted  Folding  Waveforms  after  Comparators. 


Finally,  by  combining  the  information  in  Table  2  and  Figure  11,  and  using  Equa¬ 
tion  (3.17),  we  can  obtain  a  plot  of  estimated  AOA  versus  true  AOA  as  shown  in  Figure 
12. 


Figure  12.  Estimated  AOA  versus  True  AOA. 


Next,  we  compare  the  calculated  results  with  measured  experimental  data  and 
analyze  the  differences. 
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2.  Experimental  Results  and  Analysis 

A  prototype  circuit  based  on  Figure  9  was  built  and  tested  [7,  8,  10].  The  meas¬ 
ured  AO  A  versus  true  AO  A  is  shown  in  Figure  13. 


System  Transfer  Function 


Actual  AOA  (degrees) 

|  —  Estimated  AOA  Difference 


Figure  13.  Measured  RSNS  Transfer  Function.  (From  [10].) 

Although  we  find  that  the  transfer  function  shown  in  Figure  13  exhibits  the  basic 
features  of  the  simulated  (ideal)  one  shown  in  Figure  12,  the  performance  is  not  good 
enough  to  achieve  the  DF  objective.  Possible  error  sources,  such  as  comparators  and  mi¬ 
crowave  components,  have  been  discussed  in  [7,  8].  Here,  we  analyze  the  errors  from  a 
theoretical  viewpoint.  That  is,  an  error  analysis  is  perfonned  on  each  parameter  of  the 
fundamental  Equation  (3.6)  to  show  how  errors  effect  the  performance  of  the  transfer 
function. 

First  of  all,  we  examine  the  antenna  spacing  dr  If  the  antenna  spacing  has  a  large 

fabrication  tolerance  (A  d  =  (dlheoretial  ~dactual)x\ 00%)  or  if  the  spacing  does  not  follow 

Equation  (3.11),  it  may  cause  errors  similar  to  those  shown  in  Figure  13.  The  error  simu¬ 
lation  result,  shown  in  Figure  14,  exhibits  this  phenomenon  for  the  case  of  a  spacing  tol¬ 
erance  Ad  =  ±5%. 
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Figure  14.  Transfer  Function  with  Ad  =  ±5%. 


The  second  parameter  contributing  to  estimation  errors  is  the  phase  adjuster  g. 
Two  cases  of  inaccurate  phase  adjustment  were  simulated.  Figure  15  shows  the  errors 
resulting  from  an  error  tolerance  of  Ag  =  (gtheoretiai  ~  $  actual ) x  1 00%  =  5%.  Figure  16  shows 
the  errors  caused  by  mismatched  g,  which  results  in  not  mapping  the  broadside  antenna 
response  to  the  center  of  the  dynamic  range. 
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Figure  16.  Transfer  Function  with  Mismatched  Phase  Adjustment  g. 

In  this  chapter,  we  have  reviewed  the  basic  features  of  the  RSNS  and  the  relation¬ 
ship  between  RSNS  theory  and  DF  system.  An  example  of  the  previous  RSNS-based  DF 
system  design  demonstrated  the  design  procedures.  In  the  next  chapter,  we  introduce  the 
noise  model  into  a  RSNS-based  DF  system  and  examine  the  effects  of  noise  on  the  trans¬ 
fer  function. 
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IV.  IMPROVEMENT  IN  RSNS  DF  SYSTEM  DESIGN 


In  this  chapter,  we  develop  several  basic  analytical  theories  that  were  used  in  the 
simulations  of  Chapter  V.  First,  quadrature  type  demodulation  is  introduced,  which  pro¬ 
vides  more  flexibility  in  selection  of  a  DF  processing  algorithm,  and  potential  perform¬ 
ance  improvement.  Then,  the  additive  Gaussian  noise  model  is  presented  and  integrated 
into  the  previous  RSNS  theory  of  Chapter  III.  For  perfonnance  assessment,  the  Monte 
Carlo  method  was  selected,  which  is  widely  used  in  modeling  physical  random  processes. 
Finally,  a  virtual  channel  concept  is  introduced  and  employed  to  improve  the  RSNS- 
based  DF  system  in  the  presence  of  noise.  An  overall  RSNS  system  design  block  diagram 
is  suggested  at  the  conclusion  of  this  chapter. 

A.  QUADRATURE  TYPE  DEMODULATION 

1.  Complex  Representation  of  Band-pass  Signal 

Bandpass  signals  are  defined  as  signals  whose  spectrums  are  nonzero  only  within 
a  region  of  ±  fB  Hertz  around  a  carrier  frequency  fc  (=  2 7r/(Dc).  Any  bandpass  signal  can 
be  expressed  using  complex  notation  as 

xband(P’t)  =  Re{xc(p)eJ0>cr}  (4-l) 

where  p  represents  the  position  infonnation  and  xc  ( p )  is  a  vector  phasor  that  contains 
information  on  direction,  magnitude,  and  phase.  The  phasor  can  be  expressed  as 

V  {p)  =  X1  {p)  +  jxQ  ( P )  =  \xc  (/0| e^(P)  (4-2) 

where  (f){p)  represents  the  phase  tenn  depending  on  the  position. 

Using  a  trigonometric  identity,  we  can  rewrite  Equation  (4.1)  as 

xband  {P’t)  =  xi  (p)  cos  (coct)  -  xQ  (p)  sin  (coct)  (4.3) 

where  x,  ( /;)  is  defined  as  the  in-phase  component  and  xa  ( p)  is  defined  as  the  quadra¬ 
ture  component.  With  respect  to  Equation  (4.1),  we  call  xc  (p)  the  baseband  signal  in 

complex  form  (or  complex  envelope),  and  e]0>ct  the  carrier  signal  in  complex  form.  The 
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product  of  the  two  represents  the  modulated  signal.  The  real  part  of  this  product  is  the 
transmitted  signal.  In  many  problems,  it  is  convenient  to  suppress  the  time  variation  and 
represent  the  signal  at  any  point  by  the  vector  phasor  shown  in  Equation  (4.2). 

2.  Quadrature  Type  Demodulation 

As  shown  in  Equation  (4.2),  the  x7  (p)  and  xQ[p)  terms  preserve  the  phasor’s 

complete  infonnation.  Therefore,  quadrature  type  demodulation  is  widely  used  in  many 
applications  that  need  complete  signal  infonnation.  In  order  to  further  process  antenna 
received  waveforms,  we  employ  the  quadrature  type  demodulation  structure  shown  in 
Figure  17  to  replace  the  conventional  structure  (Figure  3). 

Assume  that  the  signal  received  by  antenna  n  is  expressed  as 

v  ( Pn  >  0  =  Re  { A  ( Pn ) eMt) }  (4-4) 

where  n  =  1,2,3  is  the  antenna  element  index  and  A  ( pn )  is 

A (p„ )  =  | A(pn )| eMp"]  =  4  (pn )  +  jAQ  (pn ) .  (4.5) 

Figure  17  depicts  a  two-channel  interferometer  configuration  using  quadrature  type  de¬ 
modulation. 


Figure  17.  Quadrature  Type  Demodulation  for  a  Two-Channel  Interferometer. 
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A  quadrature  local  oscillator  signal  (LOq)  is  obtained  by  shifting  the  local  oscilla¬ 
tor  signal,  whose  frequency  is  equal  to  the  carrier  frequency  coc ,  by  90  degrees  as  shown 
in  Figure  17.  That  is 

vWq  =cos(<wct+;r/2)  =  -sin(ry<7).  (4.6) 

The  received  signal  is  fed  to  the  mixer  and  mixed  with  the  LO  signal  to  get  the  in-phase 
term 


Vl(Pn’t)  =  v{p„,t)cOs((Oct) 

=  ^4{P*)  +  ^[Ai{P*)™s{2G>ct)-AQ(p„)sm(2a)ct)] 
and  mixed  with  the  quadrature  signal  to  get  the  quadrature  tenn 

vQ{pnd)  =  -v{pn,t)^{(oct) 

=  ~[Ai{Pn)s™(2cQct)  +  AQ(pn)cos(2coct)^+^AQ(pn). 
After  filtering  with  a  low-pass  filter  (LPF),  we  obtain  the  baseband  terms  which  are 


{Pn)  =  \A{Pn)=^\A{pJ[CO&(^{Pn% 

(4.9) 

{Pn)  =  \AQ{Pn)  =  ^\A{Pn)\^((l>{Pn))- 

(4.10) 

For  a  plane  wave  | A(pn  )|  =  1  and  therefore,  we  nonnalized  Equations  (4.9)  and  (4.10)  to 
obtain  the  expressions 


(Pn  )  =  COS  (</>  (pn)). 

(4.11) 

>(p„)  =  sin  ((pip,,))- 

(4.12) 

Equations  (4.1 1)  and  (4.12)  are  referred  to  as  the  in-phase  and  quadrature  terms  of  the 
signal  received  by  antenna  element  n. 

Referring  to  Figure  17,  the  phase  differential  between  antenna  elements  1  and  3, 
A$3,  is  expressed  as 


a$3  =  0(pl)-0{p3)  =  kd1sin(0) 


(4.13) 
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where  9  is  the  incident  angle  measured  from  the  perpendicular  to  the  baseline  (broad¬ 
side),  and  k  is  the  wavenumber.  The  phase  of  antenna  element  3  is  established  as  the  ref¬ 
erence  phase,  and  dx  is  the  distance  between  element  1  and  the  reference  element  3.  A 
similar  procedure  for  antenna  elements  2  and  3  produces 

A^23  =  <t>  ( p2 )  -  (/>  ( Pi )  =  kd2  sin  ( 9 ) .  (4.14) 

At  this  point,  an  in-phase  output  for  a  channel  composed  of  antenna  element  n 
and  reference  antenna  element  r  is  defined  as 

VIn  (9)  =  cos(kdn  sin(6>))  (4.15) 

where  dn  is  the  distance  between  element  n  and  the  reference  element  r.  Correspond¬ 
ingly,  a  quadrature  output  is  defined  as 

VQn(0)  =  sin(kd„  sin(0)).  (4.16) 

Notice  that  Equations  (4.15)  and  (4.16)  can  be  related  to  Equations  (4.1 1)  and  (4.12)  as 
follows: 

K  (#)  =  cos  (kdn  sin  (0))  =  cos  ( A(pnr ) 

=  cos  [(/) ( pn )) cos (pr ))  +  sin [(/) ( pn  ))sin(^(/;r))  (4.17) 

=  VInVIr+VQnVQr 

and 


vQn  i6)  =  sin (kdn  sin  (0))  =  sin  (A</>nr ) 

=  sin  {(j){p,,))  cos  (/;,.))-  cos  {(/){p„  ))sin  (</>(pr))  (4.18) 


In  a  practical  application,  both  tenns  for  a  channel  can  be  computed  in  a  channel  proces¬ 
sor  as  shown  in  Figure  18.  After  obtaining  I  and  Q  for  each  channel,  the  processor  output 
goes  through  the  same  processing  as  the  conventional  configuration  does  (Figure  3). 

Referring  to  Equation  (3.16),  a  phase  adjustment  gt  is  added  to  channel  i  in  order 
to  map  the  channel  output  to  the  center  bin  of  the  system  dynamic  range  when  9  =  0°. 
That  is, 
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(4.19) 


Vj  (0)  =  cos  ( kdj  sin [0)  +  g. ) . 

Using  a  trigonometric  identity,  we  can  rewrite  Equation  (4.19)  as  follows: 

V/  (0)  =  cos  [kdj  sin  (0)  +  g. ) 

=  cos [kdi  sin(^))cosf;.  -sin(M(.  sin(6,))sin^(  (4.20) 

=  Vtt  (0)  cos  gt  -  VQi  [0)  sin  gr 

Equation  (4.20)  shows  that  phase  adjustments  can  be  implemented  in  digital  fonn  instead 
of  using  a  hardware  phase  shifter.  That  is,  we  can  acquire  the  phase  adjustment  for  each 
channel  modulus  in  advance  and  save  it  in  memory.  Figure  18  shows  the  implementation 
based  on  Equation  (4.20). 


v1{e) 

Figure  18.  Quadrature  Type  Demodulation  for  the  Second  Channel  Including  a  Digi¬ 
tal  Phase  Adjustment. 
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B. 


NOISE  MODEL 


1.  Additive  Gaussian  Noise  Model 

The  additive  Gaussian  noise  assumption  has  been  widely  adopted  for  modeling 
the  underlying  random  noise  in  many  different  applications,  since  it  usually  reduces  the 
complexity  of  the  problem  from  both  theoretical  and  empirical  standpoints.  In  this  section, 
we  add  a  bandpass  random  process  in  the  processing  algorithm  to  model  the  noise  envi¬ 
ronment. 

Consider  a  bandpass  noise  signal  in  the  following  form: 

N(t)  =  Re{|^(t)|eMV'v}  (4.21) 

where  N(t)  is  Gaussian  noise,  |^4(t)|  is  a  random  process  with  Rayleigh  first-order  prob¬ 
ability  density  function  (PDF),  and  (p{t)  is  a  random  process  with  a  uniform  first-order 
PDF  over  the  range  ( 0,  2tt]  .  It  is  convenient  to  use  the  equivalent  form: 


N{t)  =  Nj  (t)  cos  (a>ct)  -  NQ(t)sm((oct)  (4.22) 

where  N ,  (t)  and  NQ(t )  are  given  by 


JV/(0=M(0|cosW0) 

(4.23) 

NQ(t)=  ^(0sinM0)- 

(4.24) 

Combining  Equations  (4.7),  (4.8),  (4.23),  and  (4.24),  we  obtain  a 

for  the  RSNS-based  DF  system  shown  in  Figure  17.  That  is 

Gaussian  noise  model 

V*obe{Pn’t)  =  V(Pn’t)  +  Nl(  0 

(4.25) 

VQnoise{Pn’t)  =  V{Pr,’t)  +  NQ{t) 

(4.26) 

2.  Monte  Carlo  Simulation  Method 

Monte  Carlo  (MC)  methods  are  based  on  the  use  of  probability  and  statistics  to 

model  a  random  process  by  repeated  trials.  It  was  named  by  S.  Ulam,  after  the  city  of 

Monaco,  Monte  Carlo,  where  the  primary  attractions  are  casinos  providing  games  of 

chance.  MC  methods  are  applied  in  many  cases  where  a  closed-form  statistical  analysis  is 
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mathematically  complex  or  too  difficult  to  perform.  In  order  to  adopt  the  MC  method  to 
estimate  the  perfonnance  of  a  RSNS-based  DF  system  in  the  presence  of  noise  and  inter¬ 
ference,  we  first  review  the  MC  algorithm.  The  pertinent  equations  are  derived  and  then 
implemented  in  MATLAB  programs  to  obtain  the  estimation.  All  MATLAB  programs 
are  listed  and  explained  in  the  Appendix. 

A  random  variable  X  that  is  used  to  describe  a  signal  S  distorted  by  noise  W  is 
expressed  as 

X  =  S  +  W  (4.27) 

where  W  is  a  zero  mean  Gaussian  random  variable  with  variance  <72  =  1 .  Suppose  that 
the  probability  of  X  >  m  for  a  given  S  is  expressed  as 

P(m)  =  P(X>m\S )  (4.28) 

where  m  is  a  constant.  In  order  to  estimate  P{m),  we  generate  a  sequence  of  statistically 
independent  and  identically  distributed  Gaussian  random  variables  w(n)  with  zero-mean 
and  unit-variance,  where  n  =  1,2,3,...,  A.  Then,  by  adding  a  constant  value  5  to  w(n),  we 
obtain  a  new  set  of  random  variables 


x(n)  =  s  +  w(n). 

Assume  that  a  new  random  variable  y(n )  is  defined  as 

f  0,  if  x(n)<m 

v  »  =  ,  ,  \ 

[1,  if  x(n)>m 


(4.29) 


(4.30) 


Then,  the  estimation  of  P(m)  is 

A  „=i 

For  the  precision  of  Equation  (4.3 1),  we  need  to  calculate  the  mean  and  the  vari¬ 
ance  of  P(m).  The  mean  of  P(m)  is  [14] 


P(  m) 


T, E  [  y  {n )]  =  ^  NP(m  )  =  P(m) 


(4.32) 
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where  £"[•]  denotes  the  expectation  operation.  The  variance  is 


(TP(m)=E 


=  E 


P(m)  —  P(m) 


P2  (m)\-P2  (m)  =  ^ P(m)[l-P(m )]. 


Next,  suppose  that  the  precision  of  the  estimation  is  expressed  as 

p[^P(m)- P(m)  <£^j>a 


where  e  is  the  relative  error,  defined  as 


P(m)  —  P(m ) 


P(m) 

and  a  is  the  desired  outcome  percentage.  By  Chebyshev’s  inequality, 


.P(m) 


P^P(m)-P(m)\<£)>\  Hm) 

Referring  to  Equation  (4.34),  it  gives 

a<  1- 

£ " 

Substituting  Equation  (4.33)  into  Equation  (4.37),  we  obtain 

E(/n)[:i-E(ffl)] 

(1  -a)£2 

which  is  the  required  number  of  trials  for  Equation  (4.34)  to  hold. 

Now  we  sum  up  the  MC  algorithm  as  follows: 

1 .  Set  up  initial  number  of  trials,  N,  based  on  the  desirable  relative  error,  t 
centage,  a. 

2.  Generate  a  sequence  of  random  numbers  {  v  («)}'  . 

3.  Compute  the  number  of  { y(n)Y  that  exceed  m,  denoted  by  Nm. 

4.  Estimate  the  probability  P{rn)  = 


(4.33) 


(4.34) 


(4.35) 


(4.36) 


(4.37) 


(4.38) 


',  and  per- 
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5.  Use  the  estimate  P(m )  in  Equation  (4.38)  to  check  if  Equation  (4.38)  still  holds. 

If  yes,  go  to  step  6;  otherwise,  calculate  the  corresponding  Aby  Equation  (4.38) 

then  go  back  to  step  2. 

6.  Report  P(m),  then  end  the  program. 

In  the  next  chapter,  this  algorithm  is  used  to  assess  the  RSNS-based  DF  system’s 
probability  of  correctly  estimating  the  signal’s  incident  angle  in  the  presence  of  noise. 

C.  VARIABLE  RESOLUTION  SIGNAL  PROCESSING 

1.  Equal  Spacing  Antenna  Arrangement 

An  array  antenna  is  comprised  of  a  collection  of  elements  positioned  in  space  and 
interconnected.  By  the  way  of  electronic  scanning,  array  antennas  can  achieve  better  per¬ 
formance  than  mechanical  scanning  antennas.  For  this  reason,  array  antennas  play  an  im¬ 
portant  role  in  modern  EM  applications. 

The  fundamental  configuration  of  array  antennas  is  the  linear  equispaced  array  in 
which  radiation  elements  lie  equispaced  along  a  straight  line  as  shown  in  Figure  19.  Since 
this  configuration  is  simple  to  implement  with  hardware,  the  linear  equispaced  array  an¬ 
tennas  are  more  widely  used  than  any  other  array  configuration.  Therefore,  it  would  be  a 
major  advantage  if  the  RSNS  processing  technique  could  be  applied  to  periodically 
spaced  elements.  If  this  were  possible,  then  the  RSNS  method  could  be  incorporated  into 
existing  antenna  arrays,  and  there  would  be  no  requirement  to  design  a  custom  array  with 
unequally  spaced  elements. 
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From  Equation  (3.11),  when  the  frequency  is  constant,  we  obtain 


Equally  spaced  elements  require 


,  1 

df  — . 
m . 


cl  =  idx 


(4.39) 


(4.40) 


where  i  =  2,3,4,...,  N  for  N  channels.  Combining  Equations  (4.39)  and  (4.40),  we  obtain 


m , 


■  0=  i. 


m: 


(4.41) 


Since  the  RSNS  is  based  on  a  number  of  pairwise  relatively  prime  (PRP)  moduli 
[1 1],  Equation  (4.41)  does  not  hold.  The  only  way  to  satisfy  Equation  (4.41)  is  to  extend 
the  RSNS  definition  to  non-PRP  moduli.  The  cost  is  a  decrease  in  the  dynamic  range  of 
the  RSNS.  Table  4  shows  a  comparison  of  dynamic  range  values  for  both  non-PRP  and 
PRP  RSNS. 


Non-PRP  RSNS 

PRP  RSNS 

mx 

m2 

M 

mx 

m2 

M 

3 

6 

15 

3 

7 

21 

4 

8 

19 

4 

9 

32 

5 

10 

23 

5 

9 

36 

8 

16 

35 

8 

17 

64 

Table  4.  Non-PRP  versus  PRP  RSNS 


2.  Virtual  Channel 

Consider  the  two-channel  interferometer  with  quadrature  type  demodulation 
shown  in  Figure  17.  The  moduli  pair  of  m]  and  m2  are  selected  for  a  two-channel  RSNS 

whose  dynamic  range  is  M .  As  shown  in  Equation  (4.17)  and  (4.18),  the  corresponding 
in-phase  and  quadrature  terms  for  each  channel  are 


Vn  -  cos (ykdx  sin(6>)), 

(4.42) 

VQX  =sin(Mj  sin(<9)), 

(4.43) 
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V12  =  cos(M2  sin(<9)),  (4.44) 

VQ2  =  sin(M2  sin(#)).  (4.45) 

By  properly  combining  Equations  (4.42)  to  (4.45),  we  obtain  the  following  equations: 

V'\  =  %-%  =  cos(&(rf1+cf2)sin(0))  =  cos(Aa?1/sin(0)),  (4.46) 

F/2  =  VnVI2  +  VQXVQ2  —  cos(^(<i1  -  <f2)sin(<9))  =  cos(M2  sin(<9)),  (4.47) 

F'i  =  VQlVI2+VnV02  =sm[k(dx  +  <72)sin(#))  =  sin(M1,sin(0)),  (4.48) 

Eg2  =  fgif/2  -Vi\VQ2  =  sin(A:(J1  -d2)sin(<9))  =  sm (kd2  sin(<9)).  (4.49) 

Note  that  f)'  (/'  =  1,2)  can  be  regarded  as  the  in-phase  term  of  a  virtual  channel  with  ele¬ 
ments  spaced  d'.  Similarly,  V'Qi  is  the  quadrature  term  of  virtual  channel  with  elements 
spaced  d'. 

Now  suppose  that 


m2  =  hm] . 


(4.50) 


Then  using  Equations  (3.14),  (4.46),  and  (4.47): 


dx-dx+  d2 


A  M 


2  2  N 


1  1 
—  +  — 

V«i  "by 


A  M  h  +  \  A  M'? 


2  2  N  nu 


2  2 Nm' 


(4.51) 


d'=dl-d2 


A  M 


2  2  N 


1  1 


y  n\  m2  j 


A  M  h-\  _  AM'? 


2  2  N  nu 


2  2Nm' 


(4.52) 


which  give  the  relationships  for  the  virtual  spacing  to  a  new  modular  pair  [m[,  m2 } ,  hav¬ 
ing  the  dynamic  range  M',  and  the  scale  factor  ? .  Rearranging  Equations  (4.51)  and 
(4.52)  gives 


M(h  + 1) 

m2 


z 


mx 


(4.53) 
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(4.54) 


M(h- 1)  e  M' 

— - -Z  =  — 

m2  m2 

Then,  dividing  Equation  (4.53)  by  Equation  (4.54) 


f- 


h  + 1  _  m2 
h- 1 


(4.55) 


These  relationships  allow  RSNS  processing  to  be  perfonned  for  elements  that  do 
not  meet  the  physical  spacing  requirements  set  forth  in  Chapter  III.  The  overall  modified 
processing  design  for  this  so-called  variable  resolution  approach  is  described  in  the  next 
section. 


D.  OVERALL  DESIGN 

From  a  system  level  view,  a  RDF  device  consists  of  four  essential  elements  [2],  A 
block  diagram  of  a  generalized  RDF  device  is  shown  in  Figure  20. 


Figure  20.  Functional  Blocks  for  a  Direction  Finder.  (After  [2].) 


The  first  box  is  the  antenna,  which  extracts  electromagnetic  energy  and  converts  it 
into  a  signal  containing  the  AOA  information.  The  second  box  is  the  receiver,  which 
converts,  amplifies  and  processes  the  radio-frequency  signal  to  intennediate  frequency 
(IF)  or  baseband.  The  third  box  is  a  post-receiver  processor,  which  further  processes  the 
signal  in  order  to  extract  the  AOA  information.  The  fourth  box  is  the  display  device  that 
indicates  the  AOA. 

Note  that  the  operations  in  Figure  17  and  18  are  essentially  the  first  two  boxes  of 
Figure  20.  The  RSNS  decoding  block  in  Figure  9  is  equivalent  to  the  third  box  of  Figure 
20.  It  is  at  this  point  that  we  implement  the  virtual  channel  concept  and  integrate  it  into 
RSNS  decoding  block  in  order  to  improve  performance  without  changing  the  array  struc¬ 
ture.  The  flowchart  of  the  variable  resolution  technique  is  shown  in  Figure  2 1 . 
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Referring  to  Equation  (3.15),  the  higher  mt,  the  more  quantization  levels  are  re¬ 
quired.  Thus,  the  spacing  between  levels  is  smaller  and  therefore  the  coding  operation  is 
more  susceptible  to  noise  effects.  However,  I-Q  detection  allows  for  the  possibility  of  a 
variable  resolution  processor.  Initially  we  select  low-resolution  moduli  so  that  the  direc¬ 
tion  finder  is  reliable  in  a  low-SNR  environment.  This  provides  a  rough  estimate  of  the 
AOA  0L.  That  is,  the  bin  number  of  the  AOA  is  known,  but  the  bin  width  is  large.  Next 

we  use  high-resolution  moduli  to  obtain  a  more  accurate  angle  estimate  dH  within  the 

large  bin  width  of  the  first  estimate.  We  decide  which  estimated  AOA  is  displayed,  based 
on  the  reliability  of  its  measurement.  The  second  higher  resolution  estimate  is  only  re¬ 
ported  if  its  probability  of  being  correct  exceeds  a  predefined  threshold.  If  the  fine  esti¬ 
mate  angle  is  sufficiently  reliable,  then  it  is  reported;  otherwise  the  coarse  one  is  reported. 
With  this  approach,  under  a  low-SNR  condition,  we  still  obtain  a  coarse  AOA  estimate, 
but  when  the  SNR  is  high  a  much  more  accurate  AOA  estimate  is  possible. 

In  this  chapter,  several  theoretical  aspects  of  direction  finding,  signals  and  noise, 
and  the  basic  RSNS  DF  processing  have  been  discussed.  A  new  variable  resolution  RSNS 
processor  has  been  described.  In  the  next  chapter,  simulations  based  on  the  techniques 
discussed  in  this  chapter  are  used  to  predict  the  perfonnance  of  RSNS  direction  finder 
under  various  operating  parameters. 
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V.  SIMULATION,  RESULTS,  AND  ANALYSIS 


In  this  chapter,  we  present  computer  simulations  to  predict  the  performance  of 
RSNS-based  DF  systems  in  the  presence  of  noise.  The  approach  based  on  the  variable 
resolution  signal  processing  technique  is  demonstrated  and  compared  with  the  original 
RSNS  design. 


A.  NOISE  MODEL  FOR  A  RSNS-BASED  DF  SYSTEM 


In  this  section,  we  use  the  additive  Gaussian  noise  model  presented  in  Chapter  IV 
to  theoretically  assess  the  perfonnance  of  the  interferometer  shown  in  Figure  9. 
MATLAB  programs  using  the  Monte  Carlo  algorithm  were  developed  to  accomplish  the 
assessment.  All  MATLAB  programs  are  listed  and  briefly  explained  in  Appendix. 

At  the  outset,  we  introduce  the  simulation  parameters  that  are  used  throughout  the 
whole  chapter.  The  relative  error  s  was  set  at  3%,  and  the  desired  percentage  a  was  set 
at  90%.  Therefore,  the  required  number  of  trials  was  calculated  by  Equation  (4.38)  and 
dependent  on  the  corresponding  value  of  probability.  Noise  was  added  to  the  folding 
waveforms  out  of  each  channel  processor,  which  were  in  the  range  —  1 V  <  Vn  <  +1 V  and 

-IV  <  VQi  <  +1 V,  where  /  =  1,2  is  the  channel  index.  The  variance  of  the  noise  was  de¬ 
termined  from  the  specified  SNR  by 


V 

SNR  =  — Vr 

2cr 

where  Vj  is  channel  output  defined  in  Equation  (4.20). 


(5.1) 


After  noise  was  added,  the  channel  output  voltage  was  quantized  and  the  AOA 
bin  computed.  The  estimated  AOA  bin  was  obtained  by  selecting  the  bin  with  the  highest 
probability.  The  corresponding  estimated  AOA  was  obtained  by  converting  the  estimated 
bin  to  angle  units  using  Equation  (3.17).  If  the  quantized  output  value  could  not  be  en¬ 
coded  into  any  bins,  no  estimated  AOA  bin  was  obtained,  and  the  estimated  angle  was  set 
at  zero.  The  probability  of  a  correct  bin  was  defined  as  the  probability  of  the  estimated 
bin  being  equal  to  the  theoretical  bin.  The  probability  of  a  bin  error  was  defined  as  one 

hundred  percent  minus  the  probability  of  the  correct  bin. 
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First,  we  examine  the  estimated  AOA  versus  true  AOA  for  moduli  {8,17}  with 

the  SNR  being  set  at  20  dB.  Figure  22  shows  the  simulation  results.  Observing  Figure  22, 
we  find  that  the  transfer  function  “fades”  in  two  symmetrical  regions  about  0°.  That  is, 
under  the  same  operating  condition,  some  spatial  regions  may  be  estimated  incorrectly 
with  higher  probability.  This  observation  may  be  helpful  in  developing  a  solution  to  the 
“fading”  problem. 


Figure  22.  Transfer  Function  for  Moduli{8,17},  £  =  2/V3  ,  and  SNR  =  20dB. 

Next,  we  investigated  the  probability  of  correct  bin  estimation  over  a  FOV  with 
variable  SNR.  This  investigation  allowed  us  to  determine  the  RSNS  processor  operating 
environment.  Figure  23  shows  the  simulation  results.  From  this  figure,  we  find  that  it  re¬ 
quires  approximately  a  SNR  of  30  dB  to  achieve  reliable  estimation,  which  we  define  as 
the  probability  of  choosing  the  correct  bin  being  greater  than  50%. 
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Figure  23.  Probability  of  Correct  Bin  Estimation  for  Moduli  {8, 17}  and  %  =  2/V3 

with  Variable  SNR. 

Since  the  resolution  of  a  RSNS  DF  only  depends  on  two  parameters:  the  dynamic 
range  M  and  scale  factor  E,  [7-1 1],  we  simulated  the  DF  processing  for  different  moduli 
pairs  which  have  the  same  M  and  E,  to  see  how  noise  affects  their  performance.  Figures 
24  to  26  show  the  simulation  results  for  three  moduli  pairs  with  SNR  =  30  dB  and 

#  =  2/V3. 


Figure  24.  Noise  Model  for  Moduli  {5,23}  (a)  Probability  of  Bin  Error  and  (b)  Quan¬ 
tized  Output  for  Modulus  m1  =  5. 
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Figure  25.  Noise  Model  for  Moduli{8, 17)  (a)  Probability  of  Bin  Error  and  (b)  Quan¬ 
tized  Output  for  Modulus  ml  =8. 


Figure  26.  Noise  Model  for  Moduli{l0,13)  (a)  Probability  of  Bin  Error  and  (b)  Quan¬ 
tized  Output  for  Modulus  mx  =  10. 

From  the  figures,  we  find  that  the  probability  of  a  bin  error  in  all  three  cases  is 
roughly  the  same,  but  the  fading  regions  (i.e.,  areas  of  high  probability  of  bin  error)  are 
quite  different.  The  phenomenon  is  that  the  greater  the  difference  is  between  the  two 
moduli  values,  the  closer  together  the  fading  regions  appear.  Furthermore,  we  observe 
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that  the  folding  waveform  for  the  smaller  modulus  in  each  case  is  correlated  with  the 
shape  of  the  error  curve.  When  the  smaller  modulus  increases,  this  effect  is  reduced. 
Therefore,  we  conclude  that  different  moduli  with  the  same  M,  and  SNR  have  the 
same  correct  estimation  probability  but  in  different  FOV  regions. 

Not  only  is  the  probability  of  a  correct  bin  important,  but  also  the  magnitude  of 
the  error.  Figure  27  compares  the  simulation  results  for  a  noise-free  RSNS  processor  and 
one  with  additive  Gaussian  noise  (SNR  =  20 dB) .  Comparing  Figure  27  with  the  20  dB 

curve  in  Figure  23,  it  shows  that  although  some  regions  have  low  probability  of  correct 
bin  estimation,  the  processor  still  can  get  useful  AO  A  estimations.  Hence,  this  property 
may  be  helpful  for  improving  the  design.  Referring  to  Figure  13,  Figure  27  may  explain 
why  the  experimental  results  cause  huge  errors  at  some  angles. 


Angle  of  Arrival  (degrees) 

Figure  27.  Estimated  Error  Comparison  for  Moduli  {8, 17},  ^  =  2/V3 . 

B.  VARIABLE  RESOLUTION  SIGNAL  PROCESSING 
1.  PRP  RSNS 

Let  us  start  with  the  PRP  RSNS.  We  continue  using  the  moduli  ml  =  8,  and 
m2  =  17  to  illustrate  the  virtual  channel  concept  and  show  how  the  concept  increases  the 
system  resolution.  We  shall  use  a  frequency,  fc  =  8  GHz,  scale  factor,  %  =  1,  for  which 
the  FOV  is  ±90%,  and  dynamic  range,  M  =  64. 
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By  Equation  (4.50),  we  obtain  a  modular  ratio,  h  - 17/8.  Substituting  the  ratio  h 
into  Equation  (4.55),  we  get  the  new  moduli  m[  =  9  and  in',  =  25.  The  new  dynamic  range 

is  acquired  from  Equation  (3.9),  which  is  M  =  84.  Notice  that  it  is  higher  than  the  dy¬ 
namic  range  of  original  moduli  M.  By  Equation  (4.53),  the  new  scale  factor  is 

-  150/1 19  which  implies  that  the  FOV  is  about  ±53°.  Figure  28  shows  the  transfer 
function  simulated  for  the  new  design  settings. 


Figure  28.  Transfer  Function  for  Moduli{9,25}  and  =  150/119. 

Although  we  could  redesign  the  array  for  a  new  £  equal  to  we  also  must 
change  the  antenna  element  spacings.  In  the  virtual  channel  approach,  we  not  only  in¬ 
crease  resolution  with  a  higher  dynamic  range  value,  but  also  keep  the  same  physical 
element  spacings.  Therefore,  if  the  signal  arrives  within  the  new  FOV,  the  variable  reso¬ 
lution  approach  can  be  used  to  reduce  the  estimation  errors. 

At  this  point,  we  have  shown  the  new  design  approach  in  the  noise-free  operating 
condition.  We  continued  to  investigate  it  under  the  noise-added  operating  condition.  First, 
we  examined  the  probability  of  correct  bin  estimation.  Figure  29  shows  a  comparison  of 
the  two  systems  for  a  SNR  of  30  dB.  Notice  that  the  two  systems  have  slightly  different 
fading  regions  (i.e.,  low  probability  regions)  that  can  be  used  to  compensate  for  each 
other’s  fading  region. 


44 


Figure  29.  Comparison  for  Moduli {8, 17)  and  Moduli {9, 25). 

Next,  we  integrated  the  virtual  channel  concept  into  our  incident  angle  estimation 
algorithm,  whose  flowchart  is  shown  in  Figure  21.  Figure  30  shows  a  comparison  be¬ 
tween  the  two  methods.  Clearly  the  virtual  channel  has  resulted  in  smaller  error  than  the 
original  RSNS  moduli.  The  SNR  for  this  case  is  30  dB,  which  is  considered  a  relatively 
high  SNR. 


Angle  of  Arrival  (degrees) 


Figure  30.  Comparison  Between  Original  Method  and  Virtual  Channel  Method 

We  have  shown  that  the  virtual  channel  does  improve  our  direction  finding  accu¬ 
racy  for  high  SNR.  However,  as  mentioned  before,  a  low  SNR  may  drop  the  probability 
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of  correct  bin  estimation  significantly  (shown  in  Figure  23).  That  is,  if  the  operating  SNR 
is  not  high  enough,  then  the  virtual  channel  processing  may  not  contribute  much.  Figure 
3 1  demonstrates  this  case.  For  a  20  dB  SNR,  the  virtual  RSNS  processor  has  roughly  the 
same  performance  as  the  original  RSNS  processor. 


Figure  3 1 .  Comparison  Between  Original  Method  and  Virtual  Channel  Method 

Thus,  the  original  moduli  determine  the  noise  tolerance  of  the  system.  It  is  desir¬ 
able  to  have  moduli  providing  a  robust  estimation  in  low  SNR  environments  while  hav¬ 
ing  a  high  resolution  capability. 


2.  Non-PRP  RSNS 

The  example  in  Figure  31  belongs  to  Case  III  defined  in  Chapter  III.  We  now  ex¬ 
amine  the  other  two  cases  to  see  whether  the  virtual  channel  processing  can  be  applied. 


First  we  check  Case  I,  which  is  m,  >  3  and  m,  =  m,  + 1 .  Referring  to  Equation 
(4.55),  we  obtain 


m2  _  2  m,  + 1 
m,  1 

For  the  Case  II,  which  is  m,  >  5  and  m2  =  m,  +  2,  we  obtain 


(5.2) 


m2 

m[ 


2m,  +  2  _  m,  + 1 
2  1 


(5.3) 
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Since  both  cases  result  in  a  modulus  equal  to  one,  we  could  not  process  the  signal 
further,  unlike  Case  III.  In  other  words,  the  virtual  channel  concept  cannot  be  applied. 
However,  if  we  extend  the  RSNS  definition  to  non-PRP  moduli,  then  the  virtual  channel 
processing  can  be  applied.  Two  numerical  examples  illuminate  this  fact. 

For  Case  I,  moduli  m1  =  5  and  m2  =  6  with  scale  factor  E,  =  1  is  considered.  The 

corresponding  dynamic  range  is  M  =  27,  calculated  by  Equation  (3.7).  By  Equation  (5.1), 
we  obtain 

4  =  ^  (5.4) 

m[  2  1 

where  the  corresponding  dynamic  range  is  M  =47.  Referring  to  Equation  (4.50),  we  ob¬ 
tain  the  new  scale  factor 


f 


99 

235 


<1 


(5.5) 


which  means  FOV  >  ±90  .  That  is,  not  all  RSNS  bins  lie  in  the  visible  region.  For  this 
case,  the  effective  dynamic  range  is  M,  =  21(=33-13  +  l).  Figure  32  shows  the  results. 


Figure  32.  RSNS  Code  versus  AOA  for  Moduli  {2, 22} . 
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For  Case  II,  moduli  =  5  and  m2  =  7  with  scale  factor  £  =  1  is  considered.  The 

corresponding  dynamic  range  is  M  =  28,  calculated  by  Equation  (3.7).  By  Equation  (5.1) 
we  obtain 


m2  _  12  _  6 

"Z- T~7 


(5.6) 


where  the  corresponding  dynamic  range  is  M  =  27.  Referring  to  Equation  (4.50),  we  ob- 

32 

tain  the  new  scale  factor  <z  =  —  <1,  which  means  the  FOV  >  ±90°.  Figure  33  shows  the 

45 

results. 


Figure  33.  RSNS  Code  versus  AOA  for  Moduli  {2, 12) . 

For  both  cases,  although  the  virtual  channel  concept  still  holds,  it  does  not  con¬ 
tribute  to  a  higher  resolution  AOA  estimate.  That  is,  the  new  dynamic  ranges  are  lower 
than  original  RSNS  dynamic  range.  In  next  subsection,  we  investigate  the  possibility  of 
equal  spacings  and  analyze  the  potential  benefits. 

3.  Equi-spaced  RSNS 

As  mentioned  in  Chapter  IV,  accepting  non-PRP  RSNS  moduli  allows  equal  ele¬ 
ment  spacings,  but  at  the  expense  of  decreasing  dynamic  range.  Thus,  we  try  to  integrate 
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the  virtual  channel  processing  into  an  equi-spaced  RSNS  array  and  observe  whether  we 
can  obtain  other  benefits  to  compensate  for  the  dynamic  range  decrease. 

Consider  the  following  scenario: 

•  Two  channels  (three  elements),  where  the  elements  have  equal  spacing.  The 
frequency  is  fc=  8  GHz,  and  the  SNR  =  20  dB. 


•  A  RSNS  design  based  on  moduli  m1  =  8  and  m2  =16  with  a  scale  fac¬ 
tor  £  =  1.  The  corresponding  dynamic  range  is  M  =  35. 


By  Equation  (4.47),  we  obtain  h  =  2.  Substituting  the  ratio  h  into  Equation  (4.56), 
we  get  the  ratio  of  the  new  moduli  as 


m[  _  1 

m2  3 

Table  5  shows  several  possible  modular  pairs  which  meet  the  Equation  (5.7). 


(5.7) 


m[  2 

3 

4 

5 

6 

7 

8 

9 

10  11 

m2  6 

9 

12 

15 

18 

21 

24 

27 

30  33 

M'  15 

21 

27 

33 

39 

45 

51 

57 

63  69 

in'  i 

Table  5.  Moduli  for  —j  =  - . 

m2  3 


According  to  the  variable  resolution  processing  algorithm  shown  in  Figure  21,  ini¬ 
tially  we  select  a  low  resolution  moduli  to  get  a  coarse  AOA  estimate  which  is  more  ro¬ 
bust.  Thus,  the  moduli  m[  =  3  and  m2  =  9  are  selected.  The  noise  model  performance 


comparison  between  the  original  moduli  and  the  new  moduli  is  shown  in  Figure  34.  The 
data  shows  that  the  original  moduli  have  reliable  estimations  only  over  a  small  region 
while  the  new  moduli  yield  50%  or  higher  over  most  of  the  FOV.  This  result  is  expected; 


that  is,  the  new  moduli  provide  a  better  estimate  than  original  one  under  a  low  SNR  situa¬ 
tion. 
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Figure  34.  Noise  Model  for  Moduli  {3,9}  and  {8, 16} . 

Next,  let  us  examine  the  scale  factor,  which  affects  the  FOV.  By  Equation  (4.53), 
we  obtain  =  105/1 12  which  means  FOV  >  ±90°.  Figure  35  shows  the  simulation  result. 
In  this  case  all  RSNS  bins  are  used,  but  parts  of  the  first  and  last  bins  he  outside  of  the 
visible  region. 


Figure  35.  RSNS  Code  versus  AO  A  for  Moduli  {3, 9} . 

Next  we  examine  another  case  of  interest.  That  is  virtual  moduli  that  yield  a  reso¬ 
lution  higher  than  the  original  moduli.  The  moduli  m[  =  8  and  m2  =  24  are  chosen,  with  a 
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SNR  of  20  dB.  The  noise  model  performance  comparison  between  the  original  moduli 
and  the  new  moduli  is  shown  in  Figure  36. 


Figure  36.  Noise  Model  for  Moduli  {8,16}  and  {8,24}. 

Figure  36  again  demonstrates  the  phenomenon  that  the  RSNS  processor  with  a 
higher  dynamic  range  suffers  from  noise  more  easily  than  one  with  lower  dynamic  range. 
Notice  that  for  moduli  m[  =  8  and  m'2  =  24,  the  corresponding  scale  factor  is  q  =  105/102 

which  means  FOV  =  ±74  . 

The  above  three  RSNS  designs  share  the  same  element  spacing,  but  each  one  has 
a  different  level  of  resolution  and  robustness.  We  can  employ  all  three  moduli  at  the  same 
time  or  only  select  one  of  them  depending  on  the  operating  SNR  without  changing  the 
physical  element  spacings.  The  moduli  m[  =  3  and  m'2  =  9  provide  the  most  robust  esti¬ 
mations  at  lowest  SNR.  The  moduli  m1  =  8  and  m2=  16  provide  an  intermediate  level  of 
estimation,  while  the  moduli  m'  =  8  and  in',  =  24  provide  the  most  accurate  estimation. 

With  single  element  spacing,  a  three-level  progressive  angle  strategy  can  be  em¬ 
ployed  following  the  flowchart  of  Figure  21.  First,  we  determine  the  error  versus  FOV  for 
the  RSNS  moduli  ml  =  8  and  m2  =  16,  which  is  shown  in  Figure  37,  and  use  it  as  the  per¬ 
formance  baseline. 
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Figure  37.  Estimated  Error  for  Moduli  {8, 16} . 


Then,  following  the  algorithm  of  Figure  21,  Level  one  estimation  is  conducted 
with  moduli  m[  =  3  and  m2  =  9.  The  simulation  result  is  shown  in  Figure  38.  Compared 
with  the  original  moduli  in  Figure  37,  Figure  38  shows  the  significant  improvement  be¬ 
tween  40°  and  60°  from  broadside. 


Angle  of  Arrival  (degrees) 


Figure  38.  Comparison  Between  Original  Method  and  Level  One 

After  Level  one  processing,  the  signal  continues  to  be  processed  by  a  second  set 
of  moduli,  ml  =  8  and  m2  =  16,  which  is  called  Level  two  estimation.  Figure  39  shows  the 


52 


result  of  this  step.  This  time  the  improvement  is  not  as  dramatic,  but  there  is  still  some 
improvement. 

The  last  level  of  processing  {Level  three)  is  done  with  moduli  m[  =  8  and 
m2  =  24.  Figure  40  shows  the  final  level  processing  results.  At  this  point,  only  the  central 

regions  are  improved  from  Level  two.  Therefore,  for  this  selected  scenario,  two  level 
processing  is  enough. 

When  the  operating  SNR  is  high  enough,  we  can  directly  employ  the  moduli 
mx  =  8  and  m,  =  16  at  Level  one  processing.  The  drawbacks  of  adopting  non-PRP  RSNS 

coding  are  compensated  for  by  this  attractive  variable  resolution  signal  processing  capa¬ 
bility. 


Angle  of  Arrival  (degrees) 


Figure  39.  Estimated  Error  versus  AO  A  for  (a)  Level  One  (b)  Level  Two 
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Figure  40.  Estimated  Error  versus  AOA  for  (a)  Level  Two  (b)  Level  Three. 

In  this  chapter,  we  have  discussed  the  RSNS-based  DF  system  in  the  presence  of 
the  additive  Gaussian  noise.  Simulation  results  corroborated  the  Chapter  IV  analysis  and 
clearly  demonstrated  the  decrease  in  performance.  Also,  we  have  shown  the  effective  re¬ 
duction  in  angle  of  arrival  estimation  error  that  can  be  achieved  with  variable  resolution 
signal  processing.  In  the  next  chapter,  we  summarize  the  results  of  this  study  and  intro¬ 
duce  some  possible  topics  for  future  research. 
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CONCLUSIONS  AND  FUTURE  WORK 


This  research  has  analyzed  the  RSNS-based  DF  system  from  the  viewpoint  of  an 
additive  Gaussian  noise  model.  Computer  simulations  showed  that  the  presence  of  noise 
distorts  the  signal  amplitudes  applied  to  the  RSNS  encoding  and  caused  a  decrease  in  the 
AOA  accuracy.  Also  the  performances  of  different  moduli  with  the  same  parameters  have 

been  assessed.  Simulations  concluded  that  different  moduli  with  the  same  M,%,  and  SNR 
have  the  same  correct  estimation  probabilities,  but  in  different  FOV  regions. 

The  possibility  of  employing  equal  element  spacing  in  the  RSNS  processing  algo¬ 
rithm  was  investigated  by  extending  the  RSNS  definition  from  PRP  to  non-PRP.  The 
drawback  of  adopting  non-PRP  RSNS  coding  is  a  decrease  in  dynamic  range.  However, 
it  was  compensated  by  the  development  of  a  variable  resolution  signal  processing  ap¬ 
proach.  This  approach  leads  to  the  virtual  channel  concept,  which  refers  to  changing  the 
moduli  values  without  changing  the  actual  physical  antenna  element  spacing.  This  attrac¬ 
tive  property  allows  the  RSNS  algorithm  to  be  implemented  into  existing  antenna  arrays 
and  only  requires  modifying  the  antenna  signal  processor. 

Computer  simulations  of  a  two-channel  RSNS-based  DF  interferometer  have 
shown  the  effectiveness  of  the  variable  resolution  signal  processing  algorithm.  That  is, 
the  multiple  levels  and  the  progressive  angle  strategy  can  overcome  the  low  SNR  restric¬ 
tions  with  single  element  spacing. 

Suggested  future  work  includes: 

•  Simulating  and  conducting  the  variable  analog-to-digital  function  in  order  to 
employ  the  variable  resolution  signal  processing. 

•  Implementing  the  variable  resolution  signal  processing  algorithm  into  a  firm¬ 
ware  architecture  to  achieve  the  objective  of  a  digital  RSNS-based  DF  sys¬ 
tem. 

•  Developing  a  more  sophisticated  noise  model,  such  as  mutual  coupling  inter¬ 
ference,  into  RSNS  fundamental  equations  in  order  to  give  a  more  realistic 
signal  processing  simulation. 
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APPENDIX: 


SIMULATION  CODES 


The  programs  in  this  thesis  are  written  in  the  MATLAB  software.  All  program 
files  are  listed  and  briefly  explained  as  follows: 

1 .  findDynamicRange.m  is  used  to  find  the  dynamic  range  and  the  discrete  states  vec¬ 
tors. 

2.  findPhi.m  is  used  to  find  the  phase  adjustment  of  the  modular  channel. 

3.  enCode.m  is  used  to  quantize  the  modular  folding  waveforms  into  discrete  states. 

4.  moduli.m  is  used  to  compute  the  estimated  AOA. 

5.  estCorBin.m  is  used  to  compute  the  probability  of  correct  bin  estimation. 

6.  VResolution.m  is  used  to  demonstrate  the  three-level  variable  resolution  signal 
processing  algorithm. 
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1 .  %====================================== 

2.  %=  Program  Name:  f indDynamicRnage .m 

3 .  %====================================== 

4 .  % 

5.  function  [RSNS] =  f indDynamicRange (modulus) 

6.  [mR  mC] =size (modulus) ; 

7.  m=sort (modulus) ; 

8.  cN=mC;  %  cN=channelNumber 

9.  for  i=l:cN 

10.  s(i)=i-l;  %  s=shift 

11.  end 

12.  period=l; 

13.  for  i=l:cN 

14.  period=LCM (period, m ( i )) ;  %  period=LCM (2*ChannelNumber*modulus) 

15.  end 

16.  period=2*cN*period; 

17.  % - 

18  .  % 

19.  for  i=l:cN 

20.  for  j=l+s(i)  :  cN*m (i) +s (i) +2 

21.  dM ( i , j ) =f loor (  (j  -s(i)-l)/cN  );  %  dM=discrete  Modulus 

22.  end 

23.  for  j =cN*m ( i ) +3+S ( i )  :  2*cN*m (i) +s (i) 

24.  dM ( i , j ) =f loor (  (2*cN*m ( i) +cN- j +s ( i ) ) /cN  ); 

25.  end 

26.  k=l; 

27.  index=l; 

28.  while  (index  <  1.5*period) 

29.  dM (  i , 2*k*cN*m (i) +s (i) +1  :  2* (k+1) *cN*m (i) +s (i)  )=... 

30.  dM  (  i,l  +  s(i)  :  2*cN*m (i) +s (i)  ); 

31.  index=2* (k+1) *cN*m (i) +s (i) ; 

32.  k=k+l; 

33.  end 

34 .  end 

35 .  dM=dM 1 ; 

36 .  gb ( : , [2 : cN+1] ) =dM ( : , [1 : CN] ) ; 

37.  [sgbr  sgbc] =size (gb) ; 

38.  gb ( : , 1) = (0 : 1 : sgbr-1) 1 ; 

39.  gc=gb; 

40.  t  =  clock;  %  calculate  processing  time. 

4 1 .  max=  0 ; 

42.  maxMatrix= [] ; 

43.  k= 1 ; 

44 .  while  1  %  negelect  all  zero  component 

45.  if  gc (k, [2 : cN+1] ) ==zeros (1 , cN) 

46.  k=k+l; 

47.  else 

48.  break 

49.  end 

50.  end 

51.  % - 

52.  j  =k+ 1 ; 

53.  for  i=k: period  %  calculate  dynamic  range. - 

54.  while  j<  1.5*period 

55.  dyld=j ; 

56.  flag=0; 

57.  while  dyld  >i 
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58.  if  gc ( j , [2 : cN+1] )  ==  gc (dyld-l, [2 :cN+l] ) 

59.  flag=l; 

60.  break 

6 1 .  end 

62.  dyld=dyld-l; 

63.  end 

64 .  if  flag 

65.  break 

66.  else 

67.  j  = j +1 ; 

68.  end 

69.  end 

70  .  h (i , 1) =gc (i  +  1 , 1)  ; 

71 .  h (i,  2) =gc ( j , 1)  ; 

72.  end 

7  3 .  max=  0 ; 

74  .  maxMatrix=  []  ; 

75.  for  i=k:period 

76.  dy=h (i, 2) -h (i , 1) +1 ; 

77.  i f  dy  >  max 

78.  max=dy  ; 

79  .  maxMatrix=gc (h (i , 1)  : h (i , 2 )  ,  [l:cN+l]  )  ; 

80.  end 

81.  end 

82.  timeRead=etime (clock,  t)  ; 

83 .  gd ( : , [2: cN+2] ) =maxMatrix ( : , [1 : cN+1] ) ; 

84.  gd ( : , 1) = (0 : 1 :max-l) 1 ; 

85.  gd (max+1 , 1) =max; 

86.  gd (max+1 , 2 ) =period; 

87 .  for  k=3 : cN+2 

88.  gd (max+1 , k) =m (k-2 ) ; 

89.  end 

90.  RSNS=gd; 

91.  % - End  of  Program - 
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1 .  %====================================== 

2.  %=  Program  Name:  findPhi.m 

3 .  %====================================== 

4 .  % 

5.  function  [phi] =f indPhi (binMatrix) 

6 .  [bMr  bMc] =  size (binMatrix) ; 

7.  cN  =bMc - 2 ;  %channel  Number 

8.  M=binMatrix (bMr , 1 ) ;  %dynamic  range  value 

9.  for  i=l:cN 

10.  m ( i ) =binMatrix (bMr , i+2 ) ; 

11.  end 

12.  index=ceil  (  M/2  ); 

13.  if  M  <=  index+cN 

14.  setPts=M; 

15.  setLth=M-index; 

16.  else 

17.  setPts=index+cN; 

18 .  setLth=cN; 

19.  end 

20.  for  i=l:cN 

21.  table=modulusTable (m ( i ) , cN) ; 

22.  temp=binMatrix ( index : setPts , i+2 ) ; 

23.  temp=temp'; 

24.  for  j =1 : 2 *cN*m ( i ) 

25.  comp=table (1 , j : j +setLth) ; 

26.  if  temp  ==  comp 

27.  if (  rem (M, 2 )  ) 

28.  phi (i) =table  (2, j )  ; 

29.  else 

3  0  .  phi (i)  =  (table (2 , j ) + table (2  ,  j  +1) ) /2  ; 

31.  end 

32 .  break 

33.  end 

34 .  end 

35.  end 

36.  % - End  of  Program - 
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%=  Program  Name:  encode. m 

o _ 

o - 

2  .  % 

3.  function  [bin]  =  encode (Modulus , MixerOp) 

4.  for  i=l: Modulus 

5  .  c_m ( i ) =  0 ; 

6 .  end 

7.  for  i=l: Modulus 

8  .  threshold (i) =cos ( (2*Modulus-2*i  +  l) *pi/ (2*Modulus) ) 

9 .  end 

10.  for  j=l: Modulus 

11.  if  (  MixerOp  >  threshold(j)  ) 

12.  c_m(j)  =  1; 

13 .  else 

14.  c_m(j)  =  0; 

15.  end 

16 .  end 

17.  bin=sum (c_m ( : ) ) ; 

18.  % - End  of  Program - 
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1 .  %====================================== 

2.  %=  Program  Name:  moduli. m 

3 .  %====================================== 

4 .  % 

5.  function  [eAOA]  =moduli (Modulilnf , angle , Limit ) 

6.  c = 3  e  8  ; 

7.  f =8000e6 ;  %radio  frequency. 

8.  lambda=c/f;  %lambda  ==  wavelength 

9.  k=2*pi/lambda; 

10.  rad=pi/180; 

11.  errorCode= - 1 ; 

12.  epsilo=0.03; 

13.  alpha=0.9; 

14.  mm=ModuliInf ( 1 : 2 ) ; 

15.  sf =ModuliInf (3  )  ;  %scale  factor.  (change  here) 

16.  snr=ModuliInf (4 )  ; 

17.  theta=angle; 

18.  thetaR=theta*rad; 

19.  % - 

20.  binMatrix=f indDynamicRange (mm) ; 

21.  [bMr  bMc] =size (binMatrix) ; 

22.  M=binMatrix (bMr , 1 ) ;  %dynamic  range. 

23.  cN=bMc-2;  %channel  number. 

24.  nf  =  (M/ (2*cN)  )  . /mm;  %nf  ==  number  of  folds  (12- 

ref  18 ) 

25.  d=nf .* (lambda/2) *sf;  %d  ==  element  spacing  distance ( scaled)  (13- 
ref  18 ) 

26.  xi=findPhi (binMatrix) ;  %Phase  Adjustment. 

27.  bin= (binMatrix (1 :M, 1) )' ; 

28.  % - 

29.  %  Assume  the  reference  antenna  phase  p_r=0 

30.  %  that  is,  V_I (p_r) =1  &  V_Q(p_r)=0 

31.  %  then  Eq(4.14)  V_ml (theta) =cos (k*d_i*sin (theta) ) =V_I (p_m) 

32.  %  Eq(4.15)  V_mQ (theta) =sin (k*d_i*sin (theta) ) =V_Q (p_m) 

33.  % - 

34.  % - 

35.  if  (sf  <1) 

3  6  .  Fov= [-90  90]; 

37.  else 

38.  Fov=asin (1/sf ) /rad; 

39.  Fov= [ ( - Fov+1 )  (Fov-1)]; 

40.  end 

41.  if  (Limit(2)  <=  Fov(l)  | |  Limit(l)  >=  Fov(2)  ) 

42.  eAOA ( 1 : 4 ) =  0 ; 

43.  return; 

44 .  else 

4  5  .  Limit (1) =max (Limit (1)  , Fov (1) )  ; 

4  6  .  Limit (2 ) =min (Limit (2 )  , Fov (2 ) )  ; 

47.  Limit=Limit*rad; 

48 .  end 

49.  for  i=l:cN 

50.  V_Limit (i , 1) =cos (  k*d(i) *sin (Limit (1)  )  )*cos(xi(i))  -  sin( 

k*d (i) *sin (Limit (1) )  ) *sin (xi ( i ) ) ; 

51.  V_Limit (i , 2 ) =cos (  k*d ( i ) *sin (Limit (2 ) )  )*cos(xi(i))  -  sin( 

k*d (i) *sin (Limit (2) )  ) *sin (xi ( i ) ) ; 

52  .  enV ( 1 , i ) =enCode (mm ( i )  , V_Limit ( i , 1 ) )  ; 

53  .  enV (2 , i ) =enCode (mm ( i )  , V_Limit ( i , 2 ) )  ; 
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54 .  end 

55.  for  i=l : 2 

56.  flag=l; 

57.  for  j  =1 : M 

58.  if (  enV(i,:)  ==  binMatrix ( j , 3 : 2+cN)  ) 

59.  binLimit (i) =binMatrix ( j , 1) ; 

60.  flag=0; 

61.  break; 

62.  end 

63 .  end 

64  .  if  flag 

65.  bin=errorCode ; 

66.  disp('out  of  FOV'); 

67.  return; 

68.  end 

69.  end 

70.  % - Noise  Free - 

71.  for  i=l:cN 

72.  V_I (i) =cos (k*d (i) *sin (thetaR) ) ;  %  I  term  of  channel  Noise  Free 

73.  V_Q  (i)  =sin  (k*d  (i)  *sin  (thetaR)  )  ;  %  Q  term  of  channel  Noise  Free 

74 .  end 

75.  % - 

76.  if  (sf  <  1) 

77.  sf =1 ; 

78 .  end 

79.  N  = 1 0 0 ;  Inumber  of  trials 

80.  % - Additive  Noise  - 

81.  while  (1) 

82.  for  i=l : cN 

83.  VI (i, : ) =V_I (i) *ones (1,N) ; 

84.  VQ (i, : ) =V_Q (i) *ones (1,N) ; 

85.  end 

86  .  COmVl  =  complex (VI (1,  : )  , VQ (1,  : )  )  ; 

87  .  comV2  =  complex (VI (2, : ) , VQ (2, : ) ) ; 

88.  NcomVl  =  awgn (comVl , snr) ;  %  noise  added  complex  phasor 

for  channel  1 

89.  NcomV2  =  awgn (comV2 , snr) ; 

90  .  NV (1,  : ) =cos (xi (1) ) *real (NcomVl) -sin (xi (1) ) *imag (NcomVl)  ; 

91 .  NV (2 ,  : ) =cos (xi (2 ) ) *real (NcomV2 ) -  sin (xi (2 ) ) *imag (NcomV2 )  ; 

92 .  for  n=l :N 

93 .  NenV (1 , n) =enCode (mm (1) , NV ( 1 , n) ) ; 

94  .  NenV (2 , n) =enCode (mm ( 2 )  , NV (2 , n) )  ; 

95.  end 

96  .  enCmx=NenV' ; 

97.  for  i  =  l  :N 

98.  flag=l; 

99 .  for  j  =1 :M 

100.  if (  enCmx ( i , 1 : cN)  ==  binMatrix (j , 3 : 2+cN)  ) 

101.  Nbin(i) =binMatrix(j , 1) ; 

102.  f lag=0 ; 

103.  break; 

104 .  end 

105.  end 

106.  if  flag 

107.  Nbin ( i ) =errorCode ; 

108.  end 

109.  end 
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110.  % - 

111.  % - probalility - 

112.  for  i  =  l :  M 

113.  hit ( i ) =length (  f ind (Nbin==bin (i) )  ); 

114 .  end 

115.  hit(M+l)  =length(  f ind (Nbin==errorCode)  ); 

116.  [num  posBin] =max (hit )  ; 

117.  prob=num/N; 

118.  newN=floor(  prob* (1-prob) /( (1-alpha) * (epsilo^2)  ))  ; 

119.  if  newN  >  N 

120.  N=newN; 

121 .  VI=  []  ; 

122.  VQ=  [  ]  ; 

123.  COmVl=  [  ]  ; 

124.  COmV2=  []  ; 

125.  NcomVl=  []  ; 

126.  NcomV2=  []  ; 

12  7  .  NV=  []  ; 

12  8.  NenV=  []  ; 

12  9  .  enCmx=  []  ; 

13  0.  Nbin=  []  ; 

131.  else 

132.  Mb=posBin-l; 

133.  if  (  binLimit(l)  <=  Mb  &&  Mb  <=  binLimit(2)  ) 

134.  e AOA ( 1 ) = 1 ;  %  ID 

135.  eAOA (2 ) =asin (  (2*Mb+l) / (M*sf ) -1/sf ) /rad; 

136.  binWidth= (  asin((2*Mb-M+2)/(M*sf))- 

M) / (M*  sf ) )  ) / rad ; 

137.  eAOA ( 3 ) =  eAOA (2 ) -binWidth/2 ;  % 

138.  eAOA ( 4 ) =  eAOA ( 2 ) +binWidth/2 ;  % 

139.  else 

140.  eAOA ( 1 : 4 ) =  0 ; 

141.  end 

142.  break; 

143 .  end 

144 .  end 

145.  % - End  of  Program - 


Bin 

%  Main  AOA 
asin ( (2*Mb- 

NextLim ( 1 ) 
NextLim (2 ) 
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1 .  %====================================== 

2.  %=  Program  Name:  estCorBin.m 

3 .  %====================================== 

4 .  % 

5.  function  [Prob]  =estCorBin (Modulilnf , angle) 
c=3e8 ; 

f  =  8000e6  ;  %radio  frequency. 

6.  lambda=c/f;  %lambda  ==  wavelength 

7.  k=2*pi/lambda; 

8.  rad=pi/180; 

9.  errorCode= - 1 ; 

10.  epsilo=0.03; 

11.  alpha=0.9; 

12.  mm=ModuliInf ( 1 : 2 ) ; 

13.  sf =ModuliInf (3 ) ;  %scale  factor.  (change  here) 

14 .  snr=ModuliInf (4) ; 

15.  theta=angle; 

16.  thetaR=theta*rad; 

17.  % - 

18.  binMatrix=f indDynamicRange (mm) ; 

19.  [bMr  bMc] =size (binMatrix) ; 

20.  M=binMatrix (bMr , 1 ) ;  %dynamic  range. 

21.  cN=bMc-2;  %channel  number. 

22.  nf  =  (M/ (2*cN)  )  . /mm;  %nf  ==  number  of  folds  (12- 

ref  18 ) 

23.  d=nf .* (lambda/2) *sf;  %d  ==  element  spacing  distance ( scaled)  (13- 
ref  18 ) 

24.  xi=findPhi (binMatrix) ;  %Phase  Adjustment. 

25.  bin= (binMatrix (1 :M, 1) )' ; 

26.  % - 

27.  %  Assume  the  reference  antenna  phase  p_r=0 

28.  %  that  is,  V_I (p_r) =1  &  V_Q(p_r)=0 

29.  %  then  Eq(4.14)  V_ml (theta) =cos (k*d_i*sin (theta) ) =V_I (p_m) 


30.  %  Eq(4.15)  V_mQ (theta) =sin (k*d_i*sin (theta) ) =V_Q (p_m) 

31.  % - 

32.  % - Noise  Free - 

33.  Tbin=  []  ; 

34 .  for  i=l : cN 

35.  V_I  (i) =cos (k*d (i) *sin (thetaR) )  ;  %  I  term  of  channel  Noise  Free 

36.  V_Q  (i)  =sin  (k*d  (i)  *sin  (thetaR)  )  ;  %  Q  term  of  channel  Noise  Free 

37  .  V (i) =V_I (i) *cos (xi (i) ) -V_Q (i) *sin (xi (i) )  ; 

38.  enV(i) =enCode (mm(i) , V(i) ) ; 

39.  end 

4  0 .  for  j  =1 : M 

41.  if (  enV(l:2)  ==  binMatrix ( j , 3 : 2+cN)  ) 

42.  Tbin=binMatrix ( j , 1) ; 

43.  break; 

44  .  end 

45.  end 

46.  % - 

47.  if  (sf  <  1) 

48.  S  f  =  1 ; 

49.  end 

50.  N=100 ;  Inumber  of  trials 

51.  % - Additive  Noise  - 

52 .  while  ( 1 ) 

53  .  for  i  =  l : cN 
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54.  VI (i,  : ) =V_I (i) *ones (1,N)  ; 

55.  VQ (i, : ) =V_Q (i) *ones (1,N) ; 

56.  end 

57  .  comVl  =  complex (VI (1,  : )  , VQ (1,  : ) )  ; 

58  .  comV2  =  complex (VI (2,  : ) , VQ (2,  : ) )  ; 

59.  NcomVl  =  awgn (comVl , snr) ;  %  noise  added  complex  phasor 

for  channel  1 

60.  NcomV2  =  awgn (comV2 , snr) ; 

61 .  NV (1,  : ) =cos (xi (1) ) *real (NcomVl) -sin (xi (1) ) *imag (NcomVl)  ; 

62  .  NV (2 ,  : ) =cos (xi (2 ) ) *real (NcomV2 ) -  sin (xi (2 ) ) *imag (NcomV2 )  ; 

63.  for  n=l:N 

64  .  NenV ( 1 , n) =enCode (mm ( 1 ) , NV ( 1 , n) ) ; 

65  .  NenV (2 , n) =enCode (mm ( 2 )  , NV (2 , n) )  ; 

66.  end 

67.  enCmx=NenV' ; 

68.  for  i=l :N 

69.  flag=l; 

70.  for  j  =1 : M 

71.  if (  enCmx ( i , 1 : cN)  ==  binMatrix ( j , 3 : 2+cN)  ) 

72.  Nbin (i) =binMatrix ( j , 1) ; 

73.  flag=0; 

74 .  break; 

75.  end 

76.  end 

77.  if  flag 

78.  Nbin ( i ) =errorCode ; 

7  9  .  end 

80.  end 

81.  % - 

82.  % - probalility - 

83.  hit=length(  f ind (Nbin==Tbin)  ); 

84  .  prob=hit/N; 

85.  newN=floor(  prob* (1-prob) /( (1-alpha) * (epsilo^2) )) ; 

86  .  if  newN  >  N 

87.  N=newN; 

88  .  VI=  []  ; 

89.  VQ=  [ ]  ; 

90.  COmVl=  []  ; 

91.  comV2=[]; 

92.  NcomVl  =  []  ; 

93.  NcomV2=  []  ; 

94  .  NV=  []  ; 

95.  NenV=  []  ; 

9  6 .  enCmx= [ ]  ; 

97.  Nbin=  []  ; 

98.  else 

99.  Prob=prob*10  0  ; 

100.  break; 

101.  end 

102 .  end 

103.  % - End  of  Program - 
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1 .  %====================================== 

2.  %=  Program  Name:  VResolution . m 

3 .  %====================================== 

4 .  % 

5 .  clc 

6 .  clear 

7.  fid  =  f open (' Figure3 7v0 . txt w ') ; 

8.  modu0309= [3  9  105/112  20]; 

9.  modu0816=  [8  16  1  20]; 

10.  modu0  824=  [8  24  105/102  20]; 

thetaMatrix=-70 : 1 : 70 ; 

for  i=l : length (thetaMatrix) 

11.  angle=thetaMatrix ( i ) 

12  .  Limit= [-72  72]  ; 

13.  y ( 2 : 3 ) =  0 ; 

14.  eAOA03 09=moduli (modu03 09 , angle , Limit ) ; 

15.  eAOAref  =moduli (modu0816 , angle , Limit ) ; 

16.  y (1) =angle-eAOA0309 (2) ; 

17.  yRef=angle-  eAOAref (2); 

18.  if (eAOA0309 (1) ==1) 

19.  Limit  =  eAOA03  0  9  (3:4)  ; 

20.  eAOA0816=moduli (modu0816 , angle, Limit) ; 

21.  y (2) =angle-eAOA0816 (2) ; 

22.  i f ( e AOAO  816 ( 1 ) =  =  0 ) 

23.  y ( 2 ) =y ( 1 ) ; 

24.  y ( 3 ) =y ( 1 ) ; 

25.  else 

26.  Limit=eAOA0816 (3:4) ; 

27.  eAOA0824=moduli (modu0824 , angle , Limit) ; 

28.  y (3) =angle-eAOA0824 (2) ; 

29.  if (eAOA0824 (1) ==0) 

30.  y ( 3 ) =y ( 2 ) ; 

31.  end 

32.  end 

33 .  end 
34  . 

fprintf (fid, '%d\t%5.2f\t%5.2f\t%5.2f\t%5.2f\n' , angle, y (1 : 3 ) , yRef ) 

35.  end 

36 .  f close (fid) ; 

37  .  % 

38.  fid  =  f open (' Figure37v0 . txt ') ; 

39.  a  =  f scanf (f id,  ' %d\t%f \t%f \t%f \t%f '  ,  [5  inf]  )  ; 

40.  a  =  a ' ; 

41 .  f close (fid) ; 

42  .  [Nr  Nc] =size (a)  ; 

43. x  = (  a (  1 :Nr, 1)  ) ' ; 

44 .  yl  = (  a (  1 :Nr, 2)  ) ' ; 

45.  y2  = (  a (  1 : Nr , 3 )  ) ' ; 

4  6.  y3  = (  a (  1 : Nr ,  4 )  )  '  ; 

4  7.  y4  = (  a (  1 : Nr ,  5 )  )  '  ; 

48  .  % 

49 .  figure (1) 

50.  plot (x, y4 , ' LineWidth ' , 2 ) ; 

51.  grid; 

52.  xlabel (' Angle  of  Arrival  (degrees )',' FontSize ', 14 ) ; 

53.  ylabel (' Estimative  Error  (degrees )',' FontSize ', 14 ) ; 

54 .  set (gca,  ' YTick'  ,-60:10:60,  ' FontSize 1  ,  14)  ; 
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55.  print  -fl  -tiff  -depsc  figure37v0 
56  .  % 

57.  figure  (2) 

58.  plot (x, yl ,' LineWidth ', 2 ) ; 

59.  xlabel (' Angle  of  Arrival  (degrees )',' FontSize ', 14 )  ; 

60.  ylabel (' Estimative  Error  (degrees )',' FontSize ',  14 )  ; 

61 .  set (gca,  ' YTick'  ,-15:5:15,  ' FontSize ' , 14)  ; 

62.  grid; 

63.  print  -f2  -tiff  -depsc  figure38v0 
64  .  % 

65.  figure  (3) 

66.  subplot  (2, 1, 1)  ; 

67 .  plot (x, yl , ' LineWidth ' , 2 ) ; 

68.  xlabel (' Angle  of  Arrival  (degrees )',' FontSize ', 14 )  ; 

69.  ylabel (' Estimative  Error  (degrees )',' FontSize ',  14 )  ; 

70 .  set (gca,  'YTick' ,-15:5:15,  ' FontSize ' , 14) ; 

71.  grid; 

72.  subplot (2 , 1, 2)  ; 

73.  plot (x, y2 ,' LineWidth ', 2 )  ; 

74 .  grid; 

75.  xlabel (' Angle  of  Arrival  (degrees )',' FontSize ', 14 ) ; 

76.  ylabel (' Estimative  Error  (degrees )',' FontSize ', 14 ) ; 

77.  set (gca,  'YTick' , -15:5:15,  'FontSize' ,14) ; 

78.  print  -f3  -tiff  -depsc  figure39v0 
79  .  % 

80.  figure  (4) 

81.  subplot  (2, 1, 1)  ; 

82 .  plot (x, y2 , ' LineWidth ' , 2 ) ; 

83.  xlabel (' Angle  of  Arrival  (degrees )',' FontSize ',  14 )  ; 

84.  ylabel (' Estimative  Error  (degrees )',' FontSize ', 14 ) ; 

85 .  set (gca,  'YTick' ,-15:5:15,  ' FontSize ' , 14) ; 

86.  grid; 

87.  subplot (2 , 1, 2) ; 

88.  plot (x, y3 ,' LineWidth ', 2 ) ; 

89.  grid; 

90.  xlabel (' Angle  of  Arrival  (degrees )',' FontSize ', 14 ) ; 

91.  ylabel (' Estimative  Error  (degrees )',' FontSize ', 14 ) ; 
92  .  set (gca,  'YTick'  ,-15:5:15,  ' FontSize ' , 14)  ; 

93.  print  -f4  -tiff  -depsc  figure40v0 

94.  % - End  of  Program - 
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